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PL0T3D 

AN  INTERACTIVE  GRAPHICS  CODE  FOR  THREE  DIMENSIONAL  PLOTS 


INTRODUCTION 

he  PLOT3D  graphics  computer  code  provides  a  user  with  an  interactive  and  powerful  tool  for 
quickly  analyzing  three  dimensional  functional  data.  The  code  allows  the  user  to  choose  from  a  variety 
of  options  or  parameters  in  selecting  an  appropriate  display  of  data.  The  flexibility  of  the  code  should 
be  sufficient  to  satisfy  a  wide  range  of  user  needs.  The  code  runs  in  an  interactive  mode  so  that  the 
user  can  easily  obtain  the  desired  display  of  data.  All  options  and  parameters  (such  as  hidden  line 
removal)  are  selected  from  one  of  the  two  menus  in  the  program.  p 

- - \>  '  /lS' '  'j 

The  code  plots  the  user’s  data  as  if  the  data  were  points  on  a  three-dimensional  surface  or  solid. 
The  code  then  permits  the  user  to  change  the  orientation  from  which  the  data  are  viewed,  the  scaling  of 
the  data,  the  frequency  of  the  points  that  are  plotted  and  other  parameters  in  the  display  of  the  data. 
The  code  is  designed  so  that  a  user  will  first  develop  the  display  on  a  graphics  terminal  screen  and  then 
plot  the  data  on  a  pen  plotter. 

The  code  is  written  in  FORTRAN  77  and  can  be  used  on  any  HP  1000  computer  which  has  an  HP 
graphics  terminal  and  the  GRAPHICS  1000-II  DGL  software. 

This  computer  code  was  developed  as  part  of  the  hydrodynamics  program  of  the  Fluid  Dynamics 
Branch  and  has  been  a  useful  tool  for  the  graphical  presentation  of  results. 

1.  CODE  DESCRIPTION 

The  development  of  the  PLOT3D  graphics  computer  code  has  been  based  on  two  existing  subrou¬ 
tines;  which  perform  the  coordinate  transformation  and  plotting  of  data.  Certain  modifications  were 
necessary  to  these  routines  in  order  to  conform  with  program  and  system  requirements.  One  of  these 
subroutines  performs  the  actual  plotting  of  data  and  hidden  line  removal  and  was  originally  developed 
by  Watkins  in  a  3-D  plotting  package.1 

The  PLOT3D  code  was  implemented  on  the  HP  1000  computer  system  and  the  code  is  opera¬ 
tional  under  various  operating  systems. 

The  PLOT3D  code  is  composed  of  three  modules.  The  main  program  is  named  PLT3D  and  con¬ 
trols  the  use  of  the  two  segments  PLT3T  and  PLT3P.  The  segment  PLT3T  plots  to  the  terminal  screen 
and  provides  all  the  interactive  options.  The  segment  PLT3P  is  used  to  plot  the  last  screen  plot  on  the 
pen  plotter.  These  modules  are  described  below. 

1.1.  PLT3D 

This  is  the  main  program  of  the  computer  code  and  is  designed  as  the  driver  for  the  two  graphics 
modules.  At  the  end  of  a  graphic  task,  control  is  transferred  to  this  driver  which  determines  the  next 


Manuacript  approved  June  18,  1984. 


task  required.  The  transfer  between  main  program  and  the  other  two  modules  is  achieved  by  standard 
main-to-segment  or  segment-to-segment  control  calls.  All  input  and  program  parameters  retain  their 
values  as  they  were  last  defined.  This  allows  for  a  transfer  of  a  graphics  image  from  a  graphics  terminal 
screen  to  a  pen  plotter  without  any  additional  information. 

1.2.  PLT3T 

This  segment  of  the  code  is  the  largest  in  size  since  it  contains  most  of  the  routines  for  program 
input-output  control  plus  the  routines  for  performing  the  graphics  output  on  a  graphics  terminal.  The 
names  and  description  of  the  routines  contained  in  this  segment  follow 


PLT3T 

DMENU 

DFUNCT 

PMENU 

PFUNCT 

ADISP 

GDISP 

PLPAR 

PLOTF 

PLOTS 

PERIM 

PLTIF 

DATEP 

REFAX 

BOXIN 

TITLE 

ERRORS 

VIEWP 

DEVON 

VPMAX 

DEVOF 

LETTR 

NUMBR 


MAIN  routine  of  this  program  segment.  Provides  parameter  definition 
and  input  of  data  for  plotting 

Clears  the  screen  and  gives  the  menu  for  changing  program  parameters 

Menu  with  options  used  by  DMENU 

Clears  the  screen  and  gives  the  menu  for  program  control 

Option  menu  used  by  PMENU 

Clears  alphanumeric  display 

Clears  graphics  display 

Contains  the  interactive  input  for  all  program  parameters 
Plots  the  function  y— f(x,z)  according  to  the  specified  programs 
parameters. 

Manipulates  data  for  scaling  and  angles  of  rotation;  Calls  PLOTF 
to  plot  each  line 

Draws  the  plane  of  projection  and  perimeter  around  the  plotted  data 

Puts  data  information  on  the  same  figure  as  the  plot 

Puts  the  date  on  the  plot 

Draws  the  orientation  of  the  axes  on  the  plot 

Draws  a  box  around  the  figure 

Puts  a  figure  title  in  lower  part  of  page 

Error  output  routine 

Defines  window  and  viewport  for  graphics  output.  Also  defines 

aspect  ratio  and  character  size 

Initializes  device  and  software  for  graphics  display 

Changes  current  viewport  to  the  maximum  for  alphanumeric  output 

Turns  ofT  graphics  device 

Text  output  routine 

Numerical  text  output 


1.3.  PLT3P 


This  is  the  third  segment  of  the  code  and  it  contains  only  routines  necessary  to  perform  graphics 
output  on  a  hard  copy  device  (pen  plotter). 


NAME 

Description 

PLT3P 

Main  program  same  function  as  in  PLT3T 

ADISP 

Clears  alphanumeric  display 

PLOTS 

Same  as  in  PLT3T  segment 

PLOTF 

Plots  on  pen  plotter  instead  of  the  graphics  terminal 

PERIM 

Same  as  in  PLT3T  segment 

PLTIF 

Same  as  in  PLT3T 

DATEP 

Same  as  in  PLT3T 

REFAX 

Same  as  in  PLT3T 

ma 


BOXIN  Same  as  in  PLT3T 

TITLE  Same  as  in  PLT3T 

DEVON  Initializes  graphics  device  (pen  ploter)  and  selects  pen  color  and  speed 

ERRORS  Same  as  in  PLT3T 

VIEW  Same  as  in  PLT3T 

VPMAX  Same  as  in  PLT3T 

DEVOF  Same  as  in  PLT3T 

LETTR  Same  as  in  PLT3T 

NUMBR  Same  as  in  PLT3T 

2.  CODE  OPERATION 

Depending  upon  how  the  PLOT3D  code  has  been  installed  on  the  system  (see  Section3)  the  user 
will  access  the  program  by  the  command 

RU,  PLT3D 

The  program  displays  a  banner  to  the  terminal  and  then  calls  segment  PLT3T.  PLT3T  first 
requests  a  color  value  as  shown  in  Fig.  1 .  A  value  between  1  and  5  is  appropriate  and  will  set  the  plot¬ 
ting  color  on  the  HP-2627  terminals.  On  other  terminals  (eg.  HP-2648)  the  color  parameter  changes 
the  plotting  line  style  and  a  value  of  1  is  most  appropriate. 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


XXXXX  XXX XX 
XXXXX  *****  PLOT  3-d  *****  xxxxx 
XXX. XX  xxxxx 

XXXXX  This  oroqram  uill  oiot  a  three  dimensional  XXXXX 
XXXXX  surface  on  the  Screen  or  on  the  oen  Plotter.  XXXXX 
XXXXX  XXXXX 
XXXXX  The  default  values  of  the  olottmo  oaranetersXXXXX 
XXXXX  have  been  set  for  a  fast  riot  on  the  screen  XXXXX 
XXXXX  After  this  fast  olot  the  user  can  chanoe  the  XXXXX. 
XXXXX  oarameters  for  the  olot  desired  XXXXX 
XXXXX  XXXXX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
EXTEk  value  FOR  COLOR  i 


PLT3T  then  prompts  for  the  name  of  the  file  that  contains  the  data  to  be  plotted.  In  the  current 
program  the  file  name  must  be  six  characters  or  less.  PLT3T  then  prompts  for  the  number  of  data 
points  in  the  file  in  the  x  and  z  directions.  These  prompts  are  shown  in  Fig.  2.  The  program  reads  the 
data  to  be  plotted  in  the  following  way: 


vj 


READ  (5) 
READ  (5) 
DO  100 
READ  (5) 
CONTINUE 


(XPL(I),  I  -  I,  NXPL) 
(ZPL(J),  J  -  1,  NZPL) 

I  -  1,  NXPL 
(YPL(I,J),  J  -  1,  NZPL) 


and  the  data  must  have  been  written  to  the  file  exactly  the  same  way  and  the  maximum  values  for 
NXPL  and  NZPL  are  currently  set  to  151.  The  value  of  NXPL  does  not  have  to  be  the  same  as  the 
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value  of  NZPL.  Furthermore  the  XPL  (I)  or  ZPL  (I)  do  not  have  to  be  uniformly  spaced.  PLT3T 
then  plots  the  data  at  the  user’s  terminal  using  the  program’s  default  parameters.  Before  the  plotting  of 
data  the  program  will  prompt  for  input  of  the  axis  labels  24  characters  each  and  for  two  lines  of  figure 
caption  52  characters  each.  A  sample  plot  is  shown  in  Fig.  3.  After  each  plot  to  the  terminal,  the  pro¬ 
gram  pauses  (with  a  request  that  7"  be  entered  before  displaying  the  selection  menus). 


TVS'V 

L" 


ScrCi!=>  THE  HAnE  OF  THE  SATA  -Hi.  FROM  WHICH  44*4 
ttti  i His  PREtoKAr!  WkL  REAS  THE  DATA  TO  BE  PLOTTER  %%$t 
t*  -  ILE  tiAfii  Ic:lE  =ELlOS 

n*m***4M#******.* 

444#  SPECIF i  THE  NuhPEP  uF  DATA  POINTS  FOR  THE  #44# 
4*1*  XiNXPL 1  ARRAY.  NX? I  = 1 0 i 

**«  specify  the:  nuhber  of  data  points  for  the  «m 

44*4  2 (NZPL)  ARRAY  N2?L  =101 

Pig.  2  —  File  name  input  and  number  of  points 
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When  the  user  is  ready  to  continue,  the  program  displays  the  data  menu,  with  eleven  options,  as 
shown  in  Fig.  4.  The  parameters  that  PLT3T  uses  in  plotting  the  data  may  be  changed  by  selecting 
options  1-9.  Option  10  causes  the  display  of  the  program  menu  and  option  11  causes  control  to  return 
to  the  main  program.  Options  1-11  will  be  discussed  in  turn.  For  options  1-9,  typing  in  a  "/"  in 
response  to  the  prompt  will  preserve  the  existing  value  of  the  parameter.  After  a  value  is  entered  for  a 
parameter,  the  user  may  continue  changing  more  parameter  values  until  option  10  or  11  is  choosen. 

Selecting  option  1  clears  the  screen  and  prompts  for  entry  of  plotting  frequency  parameters  is  the 
x  and  z  directions.  These  prompts  are  shown  in  Fig.  5.  The  default  values  of  IPLX  and  1PLZ  are  2  and 
every  other  point  is  plotted.  Every  tenth  point  would  be  plotted  if  IPLX  or  IPLZ  were  given  a  value  of 
10. 


Option  2  of  the  data  menu  allows  the  user  to  reset  the  values  of  NXPL  and  NZPL  —  see  Fig.  6. 
The  number  of  points  in  either  direction  that  are  plotted  can  be  reduced  starting  always  from  point  one 
—  but  not  increased,  from  numbers  specified  when  the  file  is  read  in. 


Option  3  prompts  for  new  values  of  DELTA,  THETA  and  PHI,  as  shown  in  Fig.  7.  DELTA 
controls  the  vertical  location  of  the  plot  on  the  screen  and  normally  a  value  of  0.0  is  appropriate. 
THETA  is  the  angle  (in  degrees)  that  the  figure  is  rotated  (counter-clockwise)  around  the  x-axis.  Note 
the  reference  axis  configuration  shown  in  Fig.  3.  A  positive  value  of  THETA  rotates  the  y-axis  toward 
the  z-axis  around  the  x-axis.  Similarly  a  positive  value  of  PHI  rotates  the  z-axis  toward  the  x-axis 
around  the  y-axis. 


Figure  8  shows  the  prompts  for  changing  the  data  scaling  parameters  (Option  4).  The  default 
values  of  7,  2  and  5  will  usually  fill  the  plotting  area.  Increasing  the  values  will  enlarge  the  plot  and 
decreasing  the  values  will  decrease  the  size  of  the  plot.  If  the  size  of  the  plot  is  too  large  —  which  can 
occur  with  some  scale  factors  and  orientations  —  the  figure  would  exceed  the  allowable  plotting  area 
and  the  program  omits  the  plot.  In  such  a  case,  the  scale  factors  should  be  reduced  so  that  the  figure 
will  fit  in  the  allowable  plotting  area. 
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Key  #  i  ~  Plotting  freauencv  of  data  (2) 

Key  #  2  “  Nunber  of  data  ooints 

Key  #  3  “  Rotation  anqles  of  figure  (3S.45) 

Kev  *  4  ~  Scalinq  of  fioure  (7,2,5) 

Kev  I  S  “  Vertical  lines  for  oerineter  (0) 
Kev  *  6  "  Draw  axis  orientation  (ves) 

Key  *  7  ~  Infornation  &  Title  (ves, no) 

Kev  *  8  ~  Masking  on/off  (no,  -fast  olot) 
Key  *  9  “  Size  of  figure  (one,  i) 

Kev  #  10  ~  Continue  olottinn 

Kev  #  11  “  Exit  -  to  Main  nroqran 

Select  a  function 
Default  values  in  Daoenthesis 
Tvoe  function  nunber  Kev  *  -  l 

Fig.  4  —  Data  menu  for  parameter  input 


nn*nn%nnnn%%%nnt%nn%nt%tn*n*nnnn*%n% 

¥**#  SPECIFY  THE  FREQUENCY  OF  PLOTING  ON  TrtE  X-AXIS  **** 
nn  IF  I Pi..x  -  2,  THEN  ONE  LINE  EVERY  TUO  POINTS  ON  **** 
*.,.♦*  THE  X-AXIS  IS  PLOTED.  IPLX  =2 
n%*nnnnn*nntnnnnn%nnnnnt%ntnnnu%% 
****  SPECIFY  THE  FREQUENCY  OF  PLOT  INC  ON  THE  2-AXIS  **** 
m*  IF  IPLZ  2,  THEN  ONE  uINE  EVERY  TWO  POINTS  ON  **** 
nn  THE  Z-AXIS  is  PLOTED.  IPLZ  =2 

Fig.  S  —  Input  for  plotting  frequency  parameters 


nn%nn%nnnnn%n%Mutu%Mutn%nn%nnnMn* 

«***  SPECIF)'  THE  NUMBER  OF  DAI  A  POINTS  FOR  THE  **'*» 
#.#**  X(NXPL)  ARRAY  NXPL  -iOi 

nn  SPECIFY  THE  NUMBER  OF  DATA  POINTS  FOR  THE  **** 
n%*  Z(NZPL)  ARRAY.  NZPL  =  101 

Fig.  6  —  Input  for  number  of  data  points 


nn%n%%nnnn*nn%nnnnnun%nnnn%n*n%n* 

nn  SPECIFY  TnE  RELATIVE  VERTICAL  POSITION  C?  THE.  m* 
UU  ORIGIN  FOR  THE  AXIS  IF  DELTA  -  0  0,  THE  ORIGIN**** 
nn  IS  CENTERED  ON  TrtE  lEFT  SIDE  DELTA  -0  •} 
*ni*Mnninnnnnnnnnninnnn*nnnnnnn 
nn  SPECIFY  TiE  ROTATION  ANGLE  AROUND  THE  *-4XI3  **** 

n*t  THE  ROTATION  IS  DEFINED  CCUNTER-ClOCKWTSE  AS  ♦'*** 
nn  THE  X-AXIS  IS  POINiED  TOWARDS  YOU  ■  I  HE  iA  -35 
nntnnntnnnnnnttnnnnmnnnnnnnnnn 
nn  SPECIFY  THE  ROTATION  ANGLE  AROUND  ThE  Y-AXIS  **** 
****  THE  ROTATION  IS  DEFINED  CGUNTER-ClOCXUISE  AS  **** 
nn  THE  Y-AXIS  IS  POINTED  TOWARDS  YOU  '  PHI  -4S 


Fig.  7  —  Input  for  rotation  angles  about  x-yaxes 


iMtnnttntnnMtttnntnttnnttttMttMttnuttMt 

MU  SPECIFY  THE  SCALING  IN  THE  X-BIRECT10N  **** 

M*t  XSCALO  =7.0 

******************************************************** 

****  SPECIFY  The  SCALING  IN  THE  Y-DIRECT1QN  I*** 

****  YSCALfl  =2.0 

******************************************************** 

MM  SPECIFY  The  SCALING  IN  THE  2-DIRECTION  MM 

ttM  ZSCAlO  =5.0 

Fig.  8  —  Input  for  scaling  factor 

Figure  9  shows  the  prompts  for  option  5.  Default  values  (IPERIM  «■  1,  IPLPX  —  0,  IPLPZ  — 
0)  will  draw  a  perimeter  line  on  the  x-z  plane  indicating  the  projection  of  the  surface  plotted  onto  that 
plane.  Values  for  IPLPX  and  IPLPZ  other  than  zero  (e.g.  10)  would  cause  vertical  lines  to  be  drawn 
from  the  edge  of  the  surface  to  the  projection  plane  (every  10th  data  points)  unless  IPERIM  —  0  in 
which  case  no  perimeter  line  will  be  drawn  and  the  second  and  third  prompts  of  Fig.  9  are  omitted. 


MUUUUUMUUUMUtUUtUUUtMMtUUUttMUU 
t$ti  BO  YOU  WANT  h  PERIMETER  BRAWN  ON  YOUR  PLOTTED  ttU 
UM  FIGURE  ?  TYPE  u  FOP  NO  OR  \  CGP  YES  IPERIrt  =i 
UtMtUUUMUUUHUtUUUttMUMUUtMUtUtMtt 
MM  SPECIFY  THE  FREQUENCY  FOR  DRAWING  VERTICAL  LINES*** 

*»**  ON  THE  PERIHETcR  IF  IPLPX=0  ONLY  VERTICAL  l'INF.S*** 

*«*  AT.  THE  CORNERS  Wt'.L  BE  DRAWN  ( a-AXTS)  .  IPLPX  =0 
tMMMtMtMttUMUMttMUtUUtUUUHUUUUntM 
MU  SPECIFY  THE  FREQUENCY  FOR  DRAWING  VtPTICA..  lINES*** 

****  ON  THE  REPINE  TER  IF  IrLfZ-Q  ONLY  VERTICAL  i  JNES*** 

****  AT  THE  CORNERS  WILL  BE  DRAWN.  (Z-AXIS'.  IPLPZ  =0 

Fig.  9  —  Input  for  perimeter  parameters 

Figure  10  show  the  prompt  for  the  reference  axis  display  option.  An  example  of  the  reference 
axis  is  shown  in  Fig.  3. 


UUUttttttUtUMUtUUUUUtUUUUtUUMUMUM 
ttU  DO  rOu  WAN  ’  the  REFERENCE  Axis  DRAWN  On  ThE  MM 
%Ut  FIGURE  ?  TYPE  0  FOR  NO  OR  i  r3R  YES.  I  AXIS  =1 

Fig.  10  —  Input  for  reference  axes 


As  indicated  by  the  prompts  shown  in  Fig.  1 1  (option  7)  the  block  of  information  about  the  input 
data  at  the  top  of  the  plot  and/or  the  title  block  at  the  bottom  of  the  plot  can  be  omitted  by  entering 
values  of  0.  Values  of  1  for  INFORM  ans  IITTLE  will  retain  the  information  and  title  blocks.  For  ITI- 
TLE  -  1  the  user  is  requested  to  type  in  a  figure  title  up  to  two  lines  (52  characters  per  line). 

Figure  12  shows  the  prompt  for  option  8.  The  default  value  of  IMASK  is  0  and  hidden  lines  are 
not  removed.  Masking  of  hidden  lines  is  a  very  CPU  intensive  task  and  much  more  time  is  required  to 
draw  the  plots  if  masking  is  requested.  The  increased  plotting  time  is  very  noticeable  both  when  the 
program  is  plotting  to  the  graphics  terminal  or  to  the  pen  plotter.  Usually  masking  will  be  selected  only 
at  the  last  step  of  developing  a  presentation  plot. 


n%m*nnntnnnnnnnnn*n*nnnnn**unnn 

«***  DO  YOU  WANT  ANY  INFORMATION  Or  THE  INPUT  DATA  **n 

m*  placed  on  the  bane  pace  as  the  pi.ot  ?  inform  -i 
ttnnitnttnntntttttttnnnntnttHntnunnntu 
#*»*  DO  YOU  WANT  THE  TITLE  OF  THE  FIGURE  PLACED  ON  HU 
nn  THE  SAME  PAGE  ?  TYPE  0/NG.i/tES.  ITITLE  =1 

Fig.  1 1  —  Input  for  information  and  title  of  plot 


nmrttnnnnunmnnnnnututunnnnnnni 

****  SPECIFY  IF  HIDDEN  LINE  REMOVAL  IS  REQUIRED  IF  m* 
****  IMASK  =  &  THEN  HIDDEN  LINES  ARE  PLOTTED  IP  tx*» 
#***  IMASK  =  i  THEN  HIDDEN  LINES  ARE  NUT  PLOTTED  *»** 
***%  FOR  IMASK  =  0  A  FAST  PLOT  IS  DRAWN  IMASK  -=3 

Fig.  12  —  Input  for  hidden  line  removal 


The  prompt  for  option  9  is  shown  in  Fig.  13.  The  default  value  of  SIZE  is  1.0  and  the  program 
sizes  the  figure  to  use  the  full  available  plotting  area.  If  a  value  of  2.0  is  specified,  the  plotted  figure 
will  be  reduced  in  size  one  half  in  each  direction.  Normally  SIZE  will  be  set  to  1 .0  when  plotting  to  the 
graphics  terminal  and  to  the  small  pen  plotters  such  as  the  HP-7470.  When  plotting  to  pen  plotters 
such  as  the  HP-7475  and  the  HP-9872  a  value  of  1.4  or  2.0  should  be  used  when  the  paper  does  not 
cover  the  full  area  of  the  plotter  bed. 


nnntntnnntnnnntnnunnnntnnntnnnn 
nn  SPECIFY  THE  SIZE  OF  THE  PAGE  THAT  fOU  WANT  ID  **#» 

m*  USE.  IF  SIZE  =  1  THEN  A  FULL  SIZE  PACE  OF  THE  **** 

nn  PLOT  DEVICE  IS  USED  IF  SIZE  -  2  THEN  ONLt  ONE  ***» 

nn  HALF  GF  THE  PAGE  IS  USED.  SIZE  =i  0 

Fig.  13  —  Input  for  figure  size 


Choice  10  on  the  display  menu  brings  up  the  plotting  menu  so  that  the  changed  parameters  can  be 
used.  This  menu  is  shown  in  Fig.  14.  Option  2  replots  the  current  data  with  the  revised  parameters. 
The  user  can  change  as  many  or  as  few  parameters  between  plots  as  is  appropriate.  Option  3  and  4  will 
read  more  data  from  the  file  which  is  currently  open  and,  for  each  figure,  the  program  will  prompt  for 
the  number  of  points  in  the  x  and  z  directions.  Option  3  uses  the  default  display  parameters  and  option 
4  uses  the  current  (or  revised)  display  parameters.  Option  1  from  the  plotting  menu  causes  the  pro¬ 
gram  to  prompt  for  a  new  file  name  and  number  of  points  to  be  read  from  the  file.  It  also  resets  the 
display  parameters  to  the  default  values. 

Selecting  option  5  on  the  plotting  menu  as  well  as  option  1 1  on  the  display  menu  causes  control  to 
be  returned  to  the  main  program  and  the  prompt  in  Fig.  15  to  be  displayed.  At  this  prompt,  a  response 
of  0  will  cause  the  program  to  terminate  while  a  response  of  1  will  transfer  control  to  segment  PLT3P 
to  plot  the  current  data  on  the  pen  plotter.  If  data  are  plotted  on  the  pen  plotter,  the  program  prompts 
for  pen  color  and  pen  speed  as  shown  in  Fig.  16.  The  pen  color  numbers  refer  to  the  pen  stall  numbers 
on  the  plotter.  The  plot  frame  and  information  (if  requested)  will  be  plotted  with  the  pen  in  stall  1  and 
the  data  will  be  plotted  with  the  pen  in  the  stall  number  corresponding  to  the  pen  color  selected.  The 
program  will  pause  for  the  user  to  ready  the  plotter  before  continuing  and  drawing  the  plot. 


‘■'ey  # 


’■'tot  new  tile  -  default  oaraMeters 
Also  old  tile  -  default  oaranete?*? 


*. ev  *  2  flot  sane  tile  -  new  raranerer* 

Did  set  of  data  £>o:ntc 

*f*v  t  3  5 lot  ?a«e  tile  -  default  Dara«eterf 

New  set  o*  data  ootftf- 
nev  ♦  4  '  Plot  sane  file  -  new  oaraneters 

New  set  of  data  onints. 

><.{-"  ♦  -.  "  Exit  -  to  wain  orooraw 


Selec-  a  PL01 71*»U  function 
Type  function  nunder  dev  4-2 


\ 

Fig.  14  —  Program  menu 


nm*m***t#*m********t#*»tt*»*t*****m 

♦♦♦♦4  & o  vou  want  to  aiot  the  last  !**#♦ 
4F**4  fiotre  on  the  8f?»  olotter  ?  ♦*♦♦♦ 

*****  if  i£S  tvee  i  it  NO  tvue  a  i 

Fig.  15  —  Input  for  hard  copy  generation 


l**»  SPECIFY  THE  FEf<  DOcOP  FOP  ThE  PlOTIER  ***.♦ 

n«  ;  black  .  2  blje.  3  red.  4  green  uu 

tm  rSN  CDcGfi  -1 

U**m*tmt**mm*MM***«x:«>Mi.am*mx*»mmx** 

*»**  SPEC i^t  «E  PEN  SPEED  POP  ThE  PLOTTER.  ’1-35/  *444 

««  8ECfs*H£N:’El’  VAuOi  -0s  jUAlITY  PLOTS  c.  *U* 

»«?  PEh  SPEED  -S  o 

Fig.  16  —  Input  for  pen  plotter  parameters 


At  completion,  program  control  returns  to  the  main  segement  and  the  program  prompts  for  con¬ 
tinuation.  Again,  a  response  of  0  terminates  the  program  and  a  response  of  1  returns  control  to  seg¬ 
ment  PLT3T  so  that  the  user  can  develop  additional  figures. 

3.  CODE  INSTALLATION 

The  authors  have  been  using  the  PLOT3D  program  on  HP  1000  computers  running  under  the 
RTE  A.l  and  RTE  6/VM  operating  systems.  There  are  no  known  problems  that  would  keep  the  pro¬ 
gram  from  running  equally  well  under  the  RTE  IVB  and  RTE  A  operating  systems.  The  source  code  is 
device  indpendent,  but  with  the  use  of  GRAPHICS  1000-II  DGL  the  program  is  device  dependent  at 
load  time.  A  listing  of  the  source  code  is  given  in  Appendix  A.  Figure  17  shows  the  loader  command 
file  for  loading  PLOT3D  for  use  with  HP-2648  terminals  and  a  HP-9872S  plotter.  If  different  terminals 
and/or  plotters  are  in  use  on  the  target  system,  the  loader  command  file  will  have  to  be  changed 
accordingly. 
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Fig.  17  -  LOADER 
command  file  for  HP- 
2648  and  HP-9872S 
graphics  devices 
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APPENDIX 


FTN7X , Q 
tFILES  0,1 
$EMA(BLK1) 

PROGRAM  PLT3D<  3 , 95) ,  G.A.KERAMIDAS  30^12  =  12:93 

C*W»# 

COMMON  /BLK i/MASK (9000 ) ,X(15i) ,Z(151) . Y <151 ,151 ) 
COMMON/DATA1/XMARG , XWMAX , XINF , YINF , XL INF , YHMAX, XYMAR . ASPR 
COMNON/DATA2/I PLX , IPLZ , THETA , PHI , DELTA , XSC , YSC , ZSC , S IZE 
COMMON/DAT A3/IPLP , ISURF , INFORM, ITITLE , IPERIM, I AXIS , NFIG3 . 1MASK 
COMMQN/DATA4/XMIN , XMAX ,ZMIN , ZMAX , YMIN.YMAX 
COMMON/DAT A5/NXPL,NZPL,IPLPX,IPLPZ,XSCAL0 .YSCAL0.ZSCAL0 
CQMMQN/DATA6/INSTN , I FAST , IPLOT , IPMENU , IPEN<2) 
C0MM0N/DATA7/NHEDX,HEDX(5),NHEDZ,HEDZ(S),NHFDY,HEDY(5) 
C0MM0N/DATA8/NTITLi,TITLEi(2) ,NTITL2,TITLE2(15) ,NTITL3,TITLE3(15) 
DIMENSION  ISEGK3),ISEG2(3> 

DIMENSION  LBUF  (510) 

DATA  ANO/2HNO/AYE/2HYE/ 

DATA  ISEGl/2HPL,2HT3,2HT  / 

DATA  ISEG2/2HPL , 2HT3 , 2HP  / 

CALL  LGBUF(LBUF,5iG) 

C***l* 

cum  THIS  IS  A  DRIVER  FOR  THE  3D  PIT  PROGRAMS. 

PLT3P  IS  THE  PROGRAM  FOR  THE  PEN  PLOTTER. 

C**m  PLT3T  IS  THE  PROGRAM  FOR  ThE  TERMINAL 

cm** 

IPMENU=G 
URITEd  ,1005) 

1005  FORMAT ( ''XXXXXXXXXXXXXXXXXXXX'XXXXXXXXXXXXXXXXXXYXXXXXXXXXXXXXXXXX " 


,/,  ‘XXXAX 

.//XXXXX  *****  PLOT  3-D  ***** 

’,/,  "XXxXX 

,/ , "XXXXX  This  proaran  will  nlot  a  three  dinensionai 
,/ /XXXXX  surface  on  the  Screen  or  an  the  pen  Plotter. 


XXXXX' 

XXXXX" 

XXXXX' 

XXXXX" 

XXXXX" 

XXXXX' 


C  ,/,  "XXXXX  XXXXX' 

C  ,/,  "XXXXX  The  default  values  of  the  dotting  parawetersXXXXX ' 

C  ,/, "XXXXX  have  been  set  for  a  fast  Diet  on  the  screen  XXXXX" 

C  .//XXXXX  After  this  fast  plot  the  user  can  change  the  XXXXX' 

C  .//XXXXX  oaraneters  for  the  nlot  desired.  XXXXX" 

C  ,/ /XXXXX  XXXXX’ 

C  ,/,"XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX") 
1  CONTINUE 
IMASK=G 
IPLOT -0 
11  CONTINUE 
ID-1 

CALL  SEGLD.ISESi.IERR) 

100  CONTINUE 
URITEd, 20) 

20  FORMAT ( ’ ****«t#*t*#*#*****mt«***t**ttt*#t *********  ' 

C  ,/,"#♦*#♦  Do  ;ou  want  to  dot  the  last  T#*t*" 

C  ,//♦###♦  fiaure  on  the  PEN  nlotter  7  H) 

READ (1.1900)  A NSW 
Hi  CO  FORMAT (A2) 

TF(ANSW  EO.AYE;  GO  TO  El 

11 


.=■ 


/  / 


GO  TO  303 
31  CONTINUE 
IPLOT=i 
IP  =  2 
IPEN(i)-iO 
IPEN(2)=i 

CALL  SECLD(IS£G2.IERR) 

29 0  CONTINUE 
WRITE (1,30) 

30  FORMAT (  'tttHtttmmtttmmttttttHttltmttttfttt 
C-  ,/, "$**♦*  Do  you  want  to  plot  any  aore  ♦**!$ 
C  ,/,"♦*♦##  fiaures?  “') 

READ (l't  1000 i  ANSW 
IF(ANSW  .EQ  AYE)  CO  TO  1 
300  CONTINUE 
STOP 


,v  .v 

i * .  •  i  • 
•  ••  .,* 

.V/'1.  ■ 


FTN7X  ,0 
•FILES  0,1 
•EMA/BLKi/ 

PROGRAM  PLT3T(5) ,  G.A.KERAMIDAS  31-04.25:84 
C0MH0N/8LK 1/HASK (9000 )  .XPLdSi)  .ZPLdSl)  .YPLUSl  ,151) 
COMMON/DATA1/XMARG .XHMAX ,XINF, YINF ,XLINF , YWMAX , XYMAR , ASPR 
COMMON/DAT A2X IPLX , IPL2 .THETA , PHI , DELTA , XSC , YSC . ZSC , SIZE 
C0MM0N/DATA3/IPLP .INFORM .ITITLE ,IPERIM . IAXIS , IMASK 
C0MM0N/DATA4/XMIN . XMAX ,ZHIN, ZMAX , YMIN , YNAX 
COMMON/DATAS/NXPL .NZPL.IPLPX.IPLPZ, XSCAL0 . YSCAL0 ,  ZSC AL  Q 
COMMON/DAT A6/INSTN . IFAST, IPLOT . IPMENU 
COMMON/D AT A7/NHEDX , HEDX ( 5 ) ,NHEDZ.HEDZ<5) .NHEDY.HEDY(S) 

COMMON/DAT A8/NT ITL 1 , T ITLEi < 2 ) . NT ITL2 , T ITLE2 < IS ) . NTITL3 . T I TLE3 ( 15 ) 
DIMENSION  VIEW<4) ,IFILE<3) 

DATA  IFILE/2H  ,2H  ,2H  /IDMENU/0/ 

CALL  DEVON 

IF  (IPMENU  GE.  4)  GO  TO  210 
10  WRITEd,1010) 

1010  format ( *************** ************»**********$« ***************** • 

C  ,/,*****  SPECIFY  THE  NAME  OF  THE  DATA  FILE  FROM ‘WHICH  *«**" 
C  ,/,"****  THIS  PROGRAM  WILL  READ  THE  DATA  TO  BE  PLOTTED  ***»" 
C  ,/,*****  FILE  NAME  IFILE  =  “) 

READd  ,1011)  IFILE 
OPEN  (UNIT=5  ,FILE=IFILE  ) 

REWIND  5 

1011  FORMAT (3A2) 

15  WRITEd,101S) 

i oi 5  format ( ********************************************************* * 

C  ,/."«**  SPECIFY  THE  NUMBER  OF  DATA  POINTS  FOR  THE  *1**" 

C  ,/,"****  X(NXPL)  ARRAY.  NXPL  =  ") 

READd,*)  NXPL 
20  WRITEd,1020) 

1020  FORMAT ( ********************************************************* “ 

C  ,/,’****  SPECIFY  THE  NUMBER  OF  DATA  POINTS  FOR  THE  **#*" 

C  ,/,"****  Z(NZPL)  ARRAY.  NZPL  =  ") 

READd,*)  NZPL 
IF ( IPMENU. GE. 4)  GO  TO  120 
100  IPLX=2 
IPLZ=2 
DELTA=0 . 0 
THETA=35. 0 
PHI=45. 0 
XSCAL0=7 
YSCAL0=2 
ZSCALQ=S 
XSC=XSCAL0 
YSC=YSCAL0 
ZSC-ZSCAL0 
IPLPX=0 
IPLPZ=0 
IPERIM=1 
IAXIS=i 
INFORMS 
ITITLE=i 


-v.-v'. 

*  *  « '  %  1 
’  %  \T  \ 


■fe: 


\V-'l v* 


-V  • 

-  "  •  •  .  •  I 


*.r  •  ***•*,  «  *  •*,«'**.*  * 
V  v1/  V  .•  V.\*  **•  V 
A  . 

v’s’  o  O  •  "  .  •  •  .  • 


a' v.v.v.  vV- A  ."  .•".''A'.' 

V.  •  **  *.  •  •.  *  *.*»  *  %  •  •  ’  • 
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I  MASK =0 
SI2E=i .  0 
IFAST=0 
IQFF-l 
I0N=2 
120  CONTINUE 

READ  (5)  (XPLII) .1=1 .NXPL) 

READ  (S)  <ZPUJ).J=1.NZPL) 

DO  800  1=1, NXPL 
READ  <S)  <YPL(I >J) .J=1.N2PL) 

300  CONTINUE 
ISO  XMAX  =  XPL(l) 

XMIN  --  XMAX 
ZMAX  =  ZPL(l) 

ZMIN  =  ZMAX 
YMAX  =  YPL ( 1  1 i 
YMIN  =  YMAX 
DO  200  I  =  i.NXPi. 

DO  200  J  =  1  ,NZ?l 
IF<XPL<I) .IT  XMIN)  XMIN  =  XPLil) 
.>'F(XPL(I)  .GT.XMAX)  XMAX  =  XPKI) 
IF(ZPKJ).LT.ZMIN)  ZMIN  =  ZPuJ) 
IF(ZPL(J)  GT  ZMAX)  ZMAX  =  ZPL(J) 
IF(YPL(I,J)  lT  YMIN)  YMIN  =  YPlUJ) 
IF <YPL( I . J ) . GT  YMAX)  YMAX  =  YPl(I.J) 
200  CONTINUE 

XMARG=0  2S/SIZE 
XYMAR=XMaRG 
ASPR=1 .0 
XUHAX=i 0 . /SIZE 
YWMAX=10 ./SIZE 
X0=0.0 
Y0=0 . 0 
XPMAX=0 . 0 
YPMAX=0 . 0 
XINF=XMARG 
YINF=YUMAX-XMARG 
XLINF=i. 0/SIZE 
YLINF=i ./SIZE 
210  CONTINUE 

CALL  VIEUP (XUMAX .YWMAX .SIZE . I ERROR ) 
CALL  ZIUS( 45l , 0 ,4 ,  IDUM  .v'IEw  ,  (ERROR ) 
IF(IERROR.EG.O)  GO  TO  305 
CALL  ERRORSd  . IERROR  i 0HZ IUS  4S1  ) 
GO  TO  999 

.305  IF ( INFORM. LE  O)  GO  TO  310 
Yl=0.8*YWMAX 
YPMAX-. 2*UIEW<4) 

CALL  F'LTI  F  ( X  0 ,  XUMAX ,  Y 1 ,  YWMAX ) 

310  IF(ITITLE.Lt.Q)  GO  TO  315 
Y2=fl . 08#YUMAX 

CALL  T ITLE ( X0, XUMAX, Y0,Y2  ) 

Y0-  08*VIEW(4) 

315  CONTINUE 


CALL  ZItt5<  451 ,0,4, IDUM , VIEW . TERROR > 

!► (IERROR  EQ. 0)  GO  TO  400 

CALL  ERRORSIi .IERROR .10HZIWS.  451  ) 

GO  TO  999 
400  XV0=VIEW(i> 

YV0=yiEW(3)+YG 

XVMAX=VIEW<2) 

yvmax=view(4)-ypmax 

41  0  CALL  7yiEW(XV0.XVMAX,YVQ,YVMAX> 

CALL  ADISP 

CALL  BOX IN< .  0  , .  0  , XWMAX . YWMAX  ,00,0.0) 

CALL  PLOTS (MASK ,XPL  ,ZPL . YPL  ,NXPL ,NZPl, 1 >NXPL) 

CALL  ZMCUR 
WRITEC  i  .2000) 

RFAB( 1  ,*)  ICONT 
CALL  GDISP(IOFF) 

CALL  PMENIKIDMENU) 

IPMENU=5 

IF  (IPMENL  .EO.  11)  GO  TO  999 
CALL  PMENIK  IPMENU/ 

CALL  GO ISP (ION) 

GO  TO  (10,150,15,15,500),  IPMENU 
500  CONTINUE 

P000  FORMAT  (///z////////////////////  . 

*  "  Tvoe  <  /  >  to  continue  when  readv  :  *) 

999  CONTINUE 
CALL  OEVOF 
CALL  SEGRT  (Z99) 

END 

SUBROUTINE  PMFNUt IPMENU) 

CALL  PFUNCT 

100  CALL  ZAi..PH(3fl ,30H  Select  a  PLOTTING  function  ) 

WRITE' 1,1000) 

1000  FORMAT)”  Tvoe  function  number  Kev  ♦  =  ") 

READU.*)  IPMENU 

IF < IPMENU  LT . 1  OR. IPMENU. GT  5)  GO  TO  100 

RETURN 

END 

SUBROUTINE  PFUNCT 
CALL  ADISP 
CALL  7NEWF 

CALL  ZALPH(48,48H  Kev  *  1  ~  Plot  new  file  -  default  oaraneters  ) 
CALL  ZALPH(48,48H  Also  old  file  -  default  oaraneters  ) 


CALL  ZALPH(48,48H  Kev  ♦  2  ~  Plot  sane  file  -  new  oaraneters  ) 

CALL  ZALPH(48,48H  Old  set  of  data  ooints  ) 

CALL  ZALPH(48,48H  Kev  ♦  3  ”  Plot  sane  file  -  default  oaraneters  ) 

CALL  Z4LPri(48,48H  New  set  of  data  ooints  ) 

CALL  ZALPH(48,48H  Kev  *  4  ~  Plot  sane  file  -  new  oaraneters  ) 

CALL  ZALPH(48,48H  New  set  of  data  ooints  ) 

CALL  ZALPH(48,48H  Key  #  5  ~  Fxit  -  to  nain  orooran  ) 

CALL  ZALPHU  ,  IK  ) 

RETURN 
END 

SUBROUTINE  DMENU(IDMENU) 


I-  -.: 


10  CONTINUE 
CALL  DFUNCT 

100  CALL  ZALPHOO  ,30H  Select  a  function  ) 

CALL  ZALPH(4Q,40H  Default  values  In  papenthesis 
WRITE (1,1000) 

1000  FORMAT!"  Type  function  number  Key  ♦  =  ") 

R£AD(l,t)  I0MENU 

IFdDMENU.LT.i  .OR.IDMENU.GT.il)  GO  TO  100 
IFdDMENU  GE.10)  GO  TO  999 
CALL  PLPAR(IDMENU) 

GO  TO  10 
999  CONTINUE 
RETURN 
END 

SUBROUTINE  DFUNCT 

COMMON/DAT A2/IPLX , IPtZ .THETA . PHI . DELTA . XSC , YSC , ZSC .SIZE 
COMMON/DAT A3/IPLP , INFORM . IT ITLE . IPER IM . I AX IS , IMASK 
COMMON/DAT  A5/NXP  L . NZPL , IPLPX , IPLPZ .XSCALO , YSCAL0 .ZSCALO 
CALL  ADISP 


CALL  ZALPH(45,4SH  Kev  #  1  ~  Ploti 

CALL  ZALPH(4S,4SH  Kev  ♦  2  "*  Nunbi 

CALL  ZALPH(45.4SH  Kev  *  3  ~  Rota' 

CALL  ZALPH(4S.4Srt  Kev  «  4  ~  Seal 

CALL  ZALPH<4S,4SH  Kev  ♦  S'*  Uert 
CALL  ZALPH(4S.4SH  Kev  *  h  ‘  Draw 

CALl  7ALPH(4S,4SH  Kev  ♦  7  “  Infoi 

CALL  ZALPH14S.4SH  Kev  «  8  ~  Mask 

CALL  ZALPH(4S.4SH  Kev  »  9  -  Sue 

CALL  ZALPH(4S.4Sh  Kev  i  ID  ~ 

CALL  ZALPH(45.4SH  Kev  ♦  ii  ~  Eut 

CALL  ZALPHd  .  1H 

RETURN 

END 

SUBROUTINE  ADISP 

INTEGER  ILIST (7) .  STRING(2),  IERR 
DATA  STRING  /155S0B.1SS12&/ 


1  ~  Plottina  frepuenev  of  data  (2) 

2  "*  Nunber  of  data  oomtt 

3  “  Rotation  anales  of  fiaure  (35.45) 

4  ~  Scaling  of  fiaure  (7.2,5) 

5  v  Vertical  lines  for  oenneter  (0) 
h  ~  Draw  axis  orientation  (ves) 

7  ~  Information  4  Title  ’.ves.no) 

8  ~  Masking  on/off  (no.  fast  olot) 

9  '  Sue  of  fiaure  (one.  1) 

id  "  Continue  olottmo 

1  ~  Exit  -  to  main  iroaram 


-V-VSf.n 

.'.V.V 

.  ■  ..  . V. 
»  *  ^  * 1 


CALL  ZALPH  (4, STRING) 

RETURN 

END 


SUBROUTINE  GDISP(IGRF) 

DIMENSION  ILIST (2) 

C 

GO  TO  ( 10 ,20) .IGRF 
10  ILIST ( i )=0 

CALL  Z0ESC( 1055,1 ,0 , IlIST .DUMMY , IERROR ) 
GO  TO  9999 
20  ILIST ( 1 > =1 

CALL  ZOESC< 10SS  1.0, ILIST. DUMMY, IERR. TR> 
C 

9999  RETURN 
END 


.  •*,  •*.  **.  •*. 


..•■va-:.- *.*v  ^ 
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SUBROUTINE  PLPAR(IDMENU) 

COMMON/DATAJ VXMARG .XWMAX ,XINF , YINF ,XLINF , YWNAX .XYMAR , ASPR 
COMMGN/DATA2/IPLX , IPLZ , THETA . PHI , DELTA , XSC , YSC,ZSC . SIZE 
C0MMDN/DATA3/IPLP , INFORM , ITITLE , IPERIM , I AXIS , I MASK 
COMMON/DAT A4/XHIN , XMAX . ZHIN , ZMAX , YMIN .  YMAX 
COMMON/DAT AS/NXPL , NZPL , IPLPX . IPLPZ , XSCALO , YSCALO .ZSCALQ 
GOTO  <5,15,25,40,55,65,70.80,85).,  IDMENU 
5  WRITE ( 1 , 1005 ) 

loos  FORMATr********************************************************" 

C  ,/."**#*  SPECIFY  THE  FREQUENCY  OF  PLOTING  ON  THE  X-AXIS  It**" 

C  ,/,"****  IF  IPLX  =  2.  THEN  ONE  LINE  EVERY  TWO  POINTS  ON  ****■' 

C  ./,"****  THE  X-AXIS  IS  PLOTED.  IPLX  =  *> 

READU,*)  IPLX 
10  WRITE(1 ,1010) 

io  i  c  format < 

c  ,/,”****  SPECIFY  THE  FREQUENCY  OF  PLOTING  ON  THE  Z-AXIS  *»**“ 

C  ,/,"****  IF  IPLZ  =  2,  THEN  ONE  LINE  EVERY  TWO  POINTS  ON  *»**" 

C  ,/,"****  THE  Z-AXIS  IS  PLOTED  IPLZ  =  ") 

READU,*)  IPLZ 
GO  TO  2000 
15  URITEU  ,1015) 

1015  FORMAT ( H****|*****||**x******|||**|**||*||||||**t*t*|**|||*|****" 

C  ./.*****  SPECIFY  THE  NUMBER  OF  DAI  A  POINTS  FOP  THE  *»**“ 

C  ./,"****  X(NXPL)  ARRAY.  NXPL  =  ") 

READU.*)  NXPL 
20  WRITEU  ,1026) 

iC29  format< "********************************************************" 

C  ,/,"****  SPECIFY  THE  NUMBER  OF  DATA  POINTS  FOR  THE  It**" 

C  ,/,"#***  Z(NZPL)  ARRAY  NZPL  =  ") 

READU,*)  NZPL 
GO  TO  2000 
25  WRITEU,  1025) 

1025  FORMAT! "*xt*****«***X*t***********«*******|****|***|*********y*l" 

C  ,//****  SPECIFY  THE  RELATIVE  VERTICAL  POSITION  OF  THE  *1**" 

t.  ,/.”****  ORIGIN  FOR  THE  AXIS.  IF  DELTA  -  6.3,  THE  ORIGIN****” 

C  ,/."****  IS  CENTERED  ON  THE  LEFT  SIDE  DELTA  =  ") 

READU,*)  DEI  TA 
30  WRITE (1.1030) 

1030  FORMAT <  "********************************************************" 

C  ,/.“****  SPECIFY  THE  ROTATION  ANGLE  AROUND  THE  X-AXIS  **** " 

C  ,/,"****  THE  ROTATION  IS  DEFINED  COUNTER-CLOCKWISE  AS  ***»“ 

C  ,/."****  ThE  X-AXIS  IS  POINTED  TOWARDS  YOU  !  THETA  =  "> 

READU,*)  THETA 
3S  URITEU, 1035) 

1035  FORMAT) "*****************************************1***1**********” 

C  ,/,’****  SPECIFY  THE  ROTATION  ANGLE  AROUND  THE  Y-AXIS.  *»**" 

C  ,/,"****  THE  ROTATION  IS  DEFINED  COUNTER-CLOCKWISE  AS  ****“ 

C  ,/,”*#**  THE  Y-AXIS  IS  POINTED  TOWARDS  YOU  !  PHI  =  "> 

READU,*)  PHI 
GO  TO  2000 
40  WRITEU, 1040) 

1640  FQRMAT< “tt*tl«*tt**tf *******************************************" 

C  ./,"****  SPECIFY  THE  SCALING  IN  THE  X-DIRECTION  *»**" 
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C  ,/,"****  XSCALO  =  ■) 

READd,*)  XSCALO 
45  URITEd  ,1045) 

1045  format  ( *  ********************************************************  * 

C  ,/,"****  SPECIFY  THE  SCALING  IN  THE  Y-DIRECTION.  *1**' 

C  ,/,"****  YSCALO  =  ") 

READd,*)  YSCALO 
50  URITEd ,1050) 

1050  format  ( 'tittittttttitituMttttitttmtmttmmmmi ******* * 

C  ,/,**«*  SPECIFY  THE  SCALING  IN  THE  Z-DIRECTION.  #***' 

C  ,/,*****  ZSCALO  =  “) 

READd  ,*)  ZSCALO 
GO  TO  2000 
S5  URITEd,  1055) 

1055  FORMAT! '********************************************************' 
C  ,/,"****  DO  YOU  WANT  A  PERIMETER  DRAWN  ON  YOUR  PLOTTED  **** 

C  ,/,“****  FIGURE  ?  TYPE  0  FOR  NO  OR  1  FOR  YES.  IPERIM  =  ") 

READd,*)  IPERIM 
IF < IPER IM  LE . 0  >  GO  TO  65 
60  URITEd, 1060) 

1060  FORMAT!"********************************************************' 

C  ,/,*****  SPECIFY  THE  FREQUENCY  FOR  DRAWING  VERTICAL  LINES!**’ 

C  ,/,*****  ON  THE  PERIMETER.  IF  IPLPX=0  ONLY  VERTICAL  LINES***' 

C  ,/,"****  AT  THE  CORNERS  WILL  BE  DRAWN.  (X-AXIS) ,IPIPX  =  ") 

READd,*)  IPLPX 
URITEd, 1061) 

1061  FORMAT ( ********************************************************* ' 

C  ,/,"****  SPECIFY  THE  FREQUENCY  FOR  DRAWING  VERTICAL  LINES*** 

C  ,/,"****  ON  THE  PERIMETER.  IF  IPI.PZ=0  ONLY  VERTICAL  LINES***’ 

C  ,/,"****  AT  THE  CORNERS  WILL  BE  DRAWN.  (Z-AXIS),  IPLPZ  =  ") 

READd,*)  IPLPZ 
GO  TO  2000 
65  URITEd, 1065) 

1065  FORMAT( "********************************************************' 

C  ,/,*****  DO  YOU  WANT  THE  REFERENCE  AXIS  DRAWN  ON  THE  **** 

C  ,/,*****  FIGURE  ?  TYPE  0  FOR  NO  OR  1  FOR  YES.  IAXIS  =  ”) 

READd,*)  IAXIS 
GO  TO  2000 
70  URITEd, 1070) 

1070  FORMAT! '********************************************************' 

C  ,/,"****  DO  YOU  WANT  ANY  INFORMATION  OF  THE  INPUT  DATA  **** 

C  ,/,"****  PLACED  ON  THE  SAME  PAGE  AS  THE  PLOT  7  INFORM  =  ") 

READd,*)  INFORM 
75  WRITEd  ,1075) 

1 075  FORMAT ( ********************************************************* 
C  ,/,"****  DO  YOU  WANT  THE  TITLE  OF  THE  FIGURE  PLACED  ON  **** 

C  ,/,"****  THE  SAME  PAGE  7  TYPE  0/NO.i/YES,  ITITLE  =  "> 
READd,*)  ITITLE 
GO  TO  2000 
80  URITEd, 1080) 

1080  FORMAT ( ********************************************************* ’ 
C  ,/,“****  SPECIFY  IF  HIDDEN  LINE  REMOVAL  IS  REQUIRED.  IF  **** 

C  ,/,"****  IMASK  =  0  THEN  HIDDEN  LINES  ARE  PLOTTED.  IF  ****’ 

C  ,/,"****  IMASK  =  1  THEN  HIDDEN  LINES  ARE  NOT  PLOTTED.  **** 


•>>v- 


C  FOR  IHASK  =  0  A  FAST  PLOT  IS  DRAWN.  MASK  =  ") 

READ (1.0  IHASK 
GO  TO  2000 
35  WRITE < 1 . 1890) 

i 090  forhat <  "*tm*m**t*******mt****t**t****t****xt****imm**r 

C  ,/."****  SPECIFY  THE  SIZE  OF  THE  PAGE  THAT  YOU  WANT  TO  ***** 

C  ./,”****  USE.  IF  SIZE  =  i  THEN  A  FULL  SIZE  PAGE  OF  THE  ***»“ 

C  ./.“****  PLOT  DEVICE  IS  USED  IF  SIZE  =  2  THEN  ONlY  ONE  *«**" 

C  ,/>"****  HALF  OF  THE  PAGE  IS  USED  SIZE  =  *) 

READ<1.*>  SIZE 
IF ( SIZE .  l.E  .0.0)  SIZE=i .  0 
2600  CONTINUE 

XSC=XSCAL0/3IZE 

YSC=YSCALO/SIZE 

ZSC=ZSCALO/SIZE 

RETURN 

F-YD 

SUBROUTINE  PLOTS (MASK ,X ,Z , Y ,NX ,NZ , IDIM , IDIHEN) 

LOGICAL  LPLOTX.LPlOTZ 

DIMENSION  XX ( 151 ) ,ZZ ( iSi ) , YY ( 151 ) , VERTEX! 16 ) 

EMA  MASK (9000). XC1S1  >  ,Z(151)  .YdSKiSii 
COMMON/DATA1 /XMARG , XWMAX , XINF , Y  INF , XLINF . YWMAX , XYMAR . ASPR 
C0HM0N7DATA2/IPL X . IPLZ  . THETA , PHI . DElTA . XSC , YSC ,2SC , SIZE 
COMMON/  DATA3/IPLP .  INFORM .  ITITLE .  IPFR  !M .  IAXIS  .MASK 
C0MM0N/DATA5/NXPL .NZPL . IPLPX. iPLPZ. XSCALO , YSCAlO . ZSCALO 
rOMMON/DATA4/XMIN .XMAX .ZMIN ,ZMAX . YMIN . YMAX 
DATA  SW/-i ./ 

? 10180=  01745329252 
PIR=3. 14159 

COSPH  =  COS(  PID180  *  PHI  ) 

COSTH  =  COS<  PID180  *  THETA  ) 

STNPH  =  SIN (  PID180  *  PHI  ) 

SINTH  =  SIN C  PID180  *  THETA  > 

TMP1  =  XSC  *  COSPH 

TMP2  =  ZSC  *  SINPH 

XMX  =  AMAXK  TMP1  ,  TMP2  .  TMP1+TMP2  ,  0.  > 

XMN  =  ANINK  TMPl  .  TMP2  .  TMP1+TMP2  .  0  ) 

TMPl  =  XSC  *  SINPH  *SINTH 

TMP2  =  -ZSC  *COSPH  *SINTH 

YMX  =  AMAXK  TMPl  .  TMP2  ,  TMP1+TMP2  .  0  ) 

YMN  =  AftINK  TMPl  .  TMP2  ,  7MP1+TMP2  ,  0  ) 

TMPl  =  YSC  *  COSTH 

YMX  =  AMAXK  YMX  ,  YMX  ♦  TMPl  ) 

YMN  =  AMINK  YMN  ,  VMM  +  TMPl  ) 

XL  =  XMX  -  XMN 
YL  =  YMX  -  YMN 

IF(  XL  GT  XWMAX  OR  YL  GT  YWMAX  )  GO  TO  998 

XOR  =  XMARG  -  XMN 

XMX  =  XMARG  +  XL  ♦  XMARG 

SW=-SW 

YOR  =  ABS ( YMN ) ♦ ( YWMAX-YL ) /2 . +DELT A 
LPLOTX  =  -SINPH  *  COSTH  GT  0. 

LPI.OTZ  *  COSPH  I  COSTH  .GT.  0. 

XSCAlE  *  XSC/(XMAX-XMIN) 


.'-Tv.-.' 


YL  GT  YWMAX  )  GO  TO  998 


v  v  v  V.V-'.v 


/*>'•.** V. 

-.'  O 


YSCALE  =  YSC/<YMAX-YHIN> 

ZSCALE  =  ZSC/(ZHAX-ZNIN) 

CALL  ZMOVE (  XOR  ,  0.0) 

KPLX  =  HAX0<  IPLX  ,  1  ) 

NEXTX  =  1 

DO  450  NLINE  =  i,NX 

IF(  NLINE  .NE.  NEXTX  )  GO  TO  450 

NEXTX  =  MINfll  NLINE  +  KPLX  ,  NX  ) 

NN  =  NLINE 

IF<  LPLOTX  )  NN  =  NX  +  i  -  NLINE 
DO  400  J=i  ,  NZ 
XX(J)  =  X(NN)-XMIN 
YY(J)  =  Y(NN,J)-YMIN 
ZZ(J)  =  Z( JJ-ZMIN 
400  CONTINUE 
NN  =  NLINE 

C ALL  PLOTF (  i i 1 i . XX ,  Y Y ,  ZZ , XSCALE , Y SC ALE > ZSCALE , NN , NZ  ,  • -P HI , -THETA , 

*  XOR , YOR , XMX ,  MASK ,  VERTEX) 

450  CONTINUE 

IF  < IPERIM . GE . i ) 

*C  ALL  PER IN (X , Y ,Z , XX. YY . ZZ .XSCALE , YSCALE , ZSCALE , XMIN , YMIN . ZMIN , 

1 NX , NZ . PHI . THETA , XOR , YOR , XMX . MASK . NN , 1 ) 

KPLZ  =  MAXQ(  IPLZ  .  i  ) 

NEXTZ  =  1 

DO  550  NLINE  *  i  ,  NZ 

IF<  NLINE  .NE.  NEXTZ  )  GO  TO  550 

NEXTZ  =  MIN0<  NLINE  +  KPLZ  ,  NZ  ) 

NN  =  NLINE 

IF (  LPLOTZ  )  NN  =  NZ  +  1  -  NLINE 
DO  500  1=1.  NX 

XX(I)  =  X(I)-XMIN 
YY(I)  =  Y< I ,NN)-YMIN 
ZZ(I)  =  Z(NN)-ZMIN 
500  CONTINUE 
NN  =  NLINE 

C ALL  P LOTF (  0 1 i i , XX , Y Y . ZZ , XSCALE , YSCALE . ZSCALE ,NN , NX , -P  HI , -THETA , 

*  XOR. YOR. XMX, MASK. 0.0) 

550  CONTINUE 

IF(IPERIM.GE.l) 

*C ALL  PER IM ( X , Y , Z , XX , YY , ZZ , XSCALE , YSCALE , ZSCALE , XMIN , YMIN , ZMIN , 

1 NX , NZ , P  HI , THET A , XOR , YOR . XMX . MASK ,  NN .  2 ) 

XC0=XUMAX-4.*XMARG 
YCfl=. 851 Y UMAX 
ZCO=XCO 

XLMX=2.5*XMARG 
YLMX=3  04XMARG 
ZLMX=2.5*XMARG 

IF(IAXIS  GE.l)  CALL  REFAX<XC0,YC0,ZC0,XLMX,YLMX.ZLMX) 

XP0=XMARG 
YP0=  934YHMAX 
CALL  DATEP(XP0,YP0) 

GO  TO  999 
998  CONTINUE 

CALL  LETTRi 0.2,0 .25, . 02 ,26HFIGURE  SIZE  LIMIT  EXCEEDED, 0 . ,26  ) 


999  CONTINUE 
RETURN 
END 

SUBROUTINE  PLOTF ( IVXYZ , XDATA , YDATA , ZDATA , XSC ALE , 

IYSCALE , ZSCALE , NLINE , NPNTS , PHI , THETA , XREF , 

*YREF,XLENTH, MASK, VERTEX) 

INTEGER  HIGH,  OLDHI,  OLDLOW 

COMMON/DATA2/IPLXO,IPLZO,THETAO , PHI 0,DELTA0, XSCQ, YSCQ, ZSCO , SIZE 
C0HH0N/DATA3/IPLP , INFORH , ITITLE, IPERIH, IAXIS , IMASK 
COMNON/DATA5/NXPL,NZPL,IPLPX,IPLPZ,XSCALO,YSCALO,ZSCALO 
DIMENSION  XDATA(i),  YDATA(l) ,  ZDATA(i) , VERTEX (16) 

EMA  MASK (9000) 

DATA  INIT,  JVXYZ,  SPHI,  STHETA/-1 ,-i , -i . OE30 ,-i . 0E30/ 

IF  ( NLINE. EQ  0)  GO  TO  SSO 
IF  (NLINE. NE.i)  GO  TO  20 
PIPI=100.*SG|RT(SIZE) 

NYPI=PIPI*20 
LIMITX=(XLENTH*PIPI)+0.5 
I=LIMITX+LIMITX 
DO  10  K=i,I 
MASK(K)=INIT 
10  CONTINUE 
INIT=-1 
INCI=-1 
1-^0 

20  IF ( JVXYZ . EQ . I VXYZ )GO  TO  70 
JVXYZ=IVXYZ 
INDZ=1 
INDY=1 
INDX*i 
INDV=i 

IF(JVXYZ.LT.1OO0)GO  TO  30 
INDV=2 

JVXYZ=JVXYZ-1000 
30  IF( JVXYZ.LT  1 0 0 ) GO  TO  40 
INDX=2 

JVXYZ=JVXYZ-100 
40  IF( JVXYZ. LT  10) GO  TO  50 
INDY=2 

JVXYZ=JVXYZ-10 
50  IF(JVXY2.LT.1)GG  TO  60 
INDZ=2 

60  JVXYZ=IVXYZ 

70  IF  ( PHI. EO. SPHI. AND. THETA. EQ.3THETA)  GO  TO  30 
SPHI=SIN(  0i?4532925*PHI> 

CPHI=CQS<  0174532925IPHI) 

STHETA=SIN(.0174S32925*THETA) 

CTHETA=COS( . 01745329254THETA) 

Ali=CPHI 

A13=-SPHI 

A2i=STHETA*SPHI 

A22»CTHETA 

A23=STHETA*CPH1 

SPHI*PHI 


STHET  A=T  HETA 
8G  INCI=-INCI 

IF(I.NE.O)I=NPNTS+t 
DO  530  K=i .NPNTS 
I=I+IWCI 

GO  TO  (90,105), INDX 
90  X=XDATA(i)+(I-i)*XSCAiE 
GO  TO  no 

100  X=XDATA(I)*XSCALE 
110  GO  TO  (120, 130), INDY 
120  Y=YDATA(l)+(I-i)#YSCALE 
GO  TO  140 

130  Y=YDATA(i)*YSCALE 
14(1  GO  TO  (155,160)  ,  INDZ 
150  Z=ZDATA(i)+(NLINE-l )IZSCAL c 
GO  TO  170 

160  Z=ZDATA(I)*ZSCALE 
170  XXX=A11*X+A13#Z+XREF 
XX=XXX 

IX=(XX*PIPI)+.S 
YYY=A2i*X  +  A23*Z  +  YREF 
YY=YYY+A22*Y 
IY=( YYIPIPI )+  5 
IF  ( IX . LE . 0 ) IX=1 
IF(IX  GT  LIMITX) IX=LIrtlTX 
IF  ( IY  ,L.T.5)IY=5 
IF ( I Y . GT  NYPI )IY=NYP1 
IF  (K.NE.l)GO  TO  25G 
lGW=IX+IX 
HIGH  =LOW-i 
NlO!<i=hASK  (LOW) 

MHIGH=MASK(HICH) 

IF (HHIGH-IY) 20 0,21 0,180 
180  IF(rtL0W-IY)19G .230 ,220 
190  LOCOLD=G 
GOTO  240 

200  <1ASK'HIGH)=IY 

I F  ( Mi.OU .  EO .  —1 )  MASK  ( LOi»l ) =1 Y 
210  lOCOLD=+1 
GO  TO  240 

220  HASk  <  LOU)-IY 
230  LQCOLD=-i 

240  IF ( I MASK  >  241.241.242 

241  CALL  ZrtOVE(XX.YY) 

GO  TO  243 

242  IF(LOCOLD.NE.O)  CALL  ZHOUE!FlOAT(IX)/?IPI . FLOAT ( I Y > /P IP I ) 

243  CONTINtlF 
JX=IX 
JY=IY 
IYREC=IY 

IFdNDV.rO. l)GO  TO  530 
INDEX=INCIH 
VERTEX ( INDEX >=XX 
VERTEX(INDEX+1 ;=YY 


VERTEX  < INDEX+8 ) =XXX 

VERTEX(INDEX+9)=YYY 

IF  INLINE. NE.i)  GO  TO  530 

VERTEX(i)=XX 

VERTEX (2)=YY 

VERTEX(9)=XXX 

VERTEX(iO)=YYY 

GO  TO  530 

250  IF ( INASK. GT.O)  GO  TO  255 

CALL  ZDRAW(XX,YY) 

GO  TO  530 

255  IF(IX.NE.JX)GG  TO  260 
JY=IY 
GO  TO  280 

260  YINC=FLOAT ( I Y-JY ) / ABS ( FLOAT ( IX- JX ) ) 

INCX=( IX-JX)/IABS( IX-JX) 

YJ=JY 

270  JX=JX+INCX 
YJ=YJ+YINC 
JY=YJ+.5 
LOU=JX+JX 
HIGH=LOU-l 
HLOU=MASK (LOU) 

NHIGH=NASK(HIGH) 

280  IF(NHIGH-JY)30fl ,300,290 
290  IF ( MLOW— JY )3i0 ,320 ,320 
300  LOC=+i 

IF (LQCOLD)360 ,370 .430 
310  LOC=fl 

IF(L0C0LD)34fl,35fl  ,330 
320  LOC=-l 

IF( LOCOLD)5iO ,450 ,440 

330  I F  <  MHIGH . LE . I YREF ) CALL  ZDR AU ( FLOAT  <  JX ) /P I P I , FLOAT ( MHIGH ) /P IP  I ) 
GO  TO  350 

340  IF(MLOW  GE. IYREFJCALL  ZDRAW(FLOAT(JX)/PIPI,FLOAT(MLOy)/PIPI) 
350  CONTINUE 
GO  TO  520 

360  IF (MLOU-I YREF >370,380 ,380 

370  IF(MHIGH-IYREF)400 ,390 ,390 

380  CALL  ZDR  AU ( FLOAT ( JX  > /P IP I , FLOAT ( HLOW ) /P IP  I > 

390  CALL  ZHOUE ( FLOAT (JX)/PIPI ? FLOAT ( MHIGH > /PIP I ) 

GO  TO  430 

400  IF(MHIGH.EQ.-i)GO  TO  430 
OLDHI=HIGH-2*INCX 
IF(NASK(OLBHI)-JY)420 ,420,410 
410  CALL  ZMOVE (FLOAT ( JX ) /P I P I . FLOAT ( J  Y ) /P IP  I ) 

GO  TO  430 
420  IP=JX-INCX 

ID=NASK(OLDHI) 

CALL  ZMOVE (FLOAT ( IP ) /P IP  I , FLOAT ( ID  > /P IP  1 ) 

430  MASK(HIGH)=JY 

IF ( MLOW . E0 . -1 ) MASK ( LOU ) = JY 
CALL  ZDR AU(FLOAT  <  JX ) /P IP  I , FLOAT ( JY ) /P IP  I ) 

GO  TO  520 
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440  :P(W*:GH-:Y*EF>46Q  ,460,450 

450  IF' MlOW-’  YffEE  >A?ii ,  470.440 

4nfl  CAi.l  ZDRAM  ( FuCAT  <  J  X ) /P IP  I  .FLOAT  (MHIGN  >/P  IP  I  i 

470  CAU  ZrtOVE ( F lGAT <  j  >.  j  / P I? I . F  LOAT ( MlDM ) /P IP  1 ) 

GO  TO  5l  0 

486  OL  Dl  Oi)=_06-?*I  NC  X 

I F •  M ASc  (  Qi.P  JW ) - J Y )  49 6 . 5 0 0 . 5 0 0 
443  CAi.L  Zi'OVE'FlCATiJO/PIF'I  .FlOAT ij  Y  J/'PIAI ) 

GO  ~Q  510 

>00  :?=;x- icx 
:C'=c.asmOlDlow) 

r «_!.  move ( float < ip ) /p :p i  . float (id)/pipi> 

5i  3  MASK  #  1 0« .  r  J  •' 

CA  L  ~DF  AW (FLOAT ( J  7 ) /P IP  I , FLOAT ( JY  > /F IP  I ) 

32  G  7  R  tP  ~  J  Y 

LQCOI.D*lOC 

IF(J).W*.TX)  GO  TO  27!) 

530  CONTINUE 

IF  (MASK)  531 ,531,532 
53:  CALL  ZHOVE(XX,Y») 

GO  TO  535 

532  CALL  ZHOVE  ( FLOAT  ( J) )  /P  IP  I  ,Fi  OAT  (JY)  /?  IP  I ) 

535  IFdNDV.EO  1/GO  TO  540 
IN0EX=-INCI+6 
VERTEX ( INDEX )=XX 
VERTEX ( INDEX+1 ) =YY 
VERTEX  (If4DEX+8>=XXX 
t/f  i?TEX  <  INDcX+?)= YYY 
IP  (N'LINE  NE.DGO  TO  540 
VERTEX ( 3 )=XX 
VERTEX(4)-YY 
VFSTEXdl  )=XXX 
VERTEX# 1?)=YYY 
540  1=1-1 
RETURN 
550  INIT=G 
RETURN 
END 

SUBROUTINE  ?ERIM(X ,Y,Z,XX.YY , ZZ.XSCAi  E . YSCalF .ZSCALE.XMIN . YHI.n 
*  ZrtlN.NX.NI .Pul  THETA,  'OP  .  rC>? . ?.MX .MASK  NN, HIM) 

COMMON /DA  T A3 / 1 3  L  ? ,  INFORM ,  ITT  TLF ,  IPESIM ,  I  A.«  1 3 . 1  MASK 
COMMON/DAT AS/NXPi  ,NZP.  PLP.<  JP.P2.X5CM>  IfCfiLO  ZSCAi  1 
DIMENSION  XX(1).YY(1).ZZ(1) 

DIMENSION  XV  (2> , YV(2) .Zv(?) 

5MA  MASK (9000 ) ,X( 151 ) . Z( 151 ) . Y ( 151 . 151 > 

1 1)  JR EC--i 

IF'.'IPl.Py. .GT.fi)  GO  TO  2 
IPLFX=NX-i 
IPLPZ=NZ-1 
2  GO  TO  (4.11, ,  1DIM 
4  DO  5  J=i ,NZ 
XX(J)=X(1>-XMIN 
YY( J)=0 . 0 
ZZ(J)*Z(J)-ZHIN 


5  CONTINUE 
DO  10  J=1)NZ,IPLPZ 
IDIREC=-IDIR£C 
XV(i)=XX(J) 

XV(2)=XV(1) 

YV(i)=0 . 0 
YV(2)=0.G 

IF ( IDIREC . EQ . — i )  YV(i)=Y(l,J)-YMIN 
IF< IDIREC. EG. +1)  YV<2)=Ytl,J>-YMIN 
ZV(i)=ZZ(J) 

ZV(2)=ZV( 1 ) 

10  CALL  P  LQTF ( 0 1 1 1 , XV , YV , ZV , XSCALE , Y SCALE ( ZSCALE , NN , 2 , -PHI , 

*  -~NE7A ,  XOR .  YOR ,  XNX ,  MASK ,  C .  I) ) 

r.A-.L  F  LO  TF  ( 0 1 1 1 ,  XX ,  YY ,  ZZ ,  XSCALE ,  YSC ALE ,  ZSCALE  ,NN , NZ , -P  HI , 

i  -TfFTA ,  XOR ,  ;T,S ,  XNX  .«ASK  >0.0) 

CO  TO  21 

} ;  !>D  15  1=1. NX 

XX  (!'=>'(  T)-X1IN 
YY ( T i~S  0 
Z7(IW.(N7)-ZMIN 
IS  CONTINUE 

DO  20  I=l.NX,IPi.PX 

direc=-jd:rec 

X!l<  D-XXd) 

XV<2)=XV(1) 

Y!1(t) -0  0 

y  v ;  5 )  -  o  o 

IF <  IDIREC . EG .  -1 )  YV ( 1  )= Y i  I , >:Z >-t NiTli 
IF  (IMPEL.  EG  i-ll  \V\.2)~Y(~  ,NZ/-Y>fIN 
7V(i>*ZZ(I> 

Z‘-,'2)=ZV(i.) 

20  "ALL  PLOT" ( 0 1 H . XV , YV , Zv > XSCALE , t  SCA  t , ZSCALE ,NN  2. -PH*. 

*  -TiETA  ; Of  .  YOP. .  XKX . HAST .  3  0  • 

r  All.  j  LOT?  ( 0 1  n  I XV  ,  YY ,  ZZ  .XSCALE ,  YSC  Al  E ,  ZSCALE  .NN.NX.-PH! . 

*  -Tr  £T(. ,  < OS ,  r OS ,  XH).  ,HA5K .  u .  6 ) 

CO  Tjj  3^ 

21  DO  25  J-i  Ki 

t;=m+i-.i 

XX(T;  =  uNX;->.PiN 
YY(J>=0.0 
ZZ(J)s:(JJ»-ZnlN 
2D  CONTINUE 

DO  30  J=l.NZ,IP  PZ 
IDIRF.C-IMPEC 
J.T=NZ+1-J 
XV(1)*XX(J) 

XV(2)-XV(1) 

YV(i)=0 . 9 
YV<2)-0.0 

IF(IDJREC.EQ.-l)  YV(i)=Y(NX,JJ)-YHIN 
IF (IDIREC. E0.+1)  YV(2)=Y(NX,JJ)-YrtIN 
ZL'(i)=ZZ(  J> 

ZV(2)=ZV(1) 

30  CALL  PLOTFtOUl, XV, YV,ZV, XSCAlE, YSCAiE, ZSCALE, NN. 2, -PHI 


*  -THETA. X3R.Y&R.XMX, MASK. 0.0) 

r AO.  ? L07F 1 0 1 1  i . XX . YY  .ZZ .  XSCALE .YSCALE ,  Z SC ALE .  NN  .NZ , -P  HI , 

%  -THETA  XOS.YDR.XM/.,  MASK.  0.0) 
no  “o  4t 
7i  DC  3S  !=i  »■ 

ii=nx+x-: 

XX ( I  •=xa:;-v»IN 
Y  i  il»=f. .  0 
ZZii  ;=Z(i>-ZM';N 

3s  coming: 

IT  40  1=1. MX  IPLPX 
II=NX+i-I 
IDIREC=-ID!REC 
XV(i)=XXtI,' 

XV(2)=XV<i) 

YU(1 )=D . 0 
YV(2)=YV(i) 

JRIDIREC.E9  -i)  YU(i )=Y( 'I . i > — YMTN 
IFIIDISEC.EQ.+l)  YV(2)=YiiM)-YMTN 
ZV(1)=ZZ(I) 

zv(2)=zy(t) 

40  CALL  PLOTFIOiii, XV, YV,ZV, XSCALE, YSCALE.7SCAIE.NK  .2. -Shi 

*  -THETA. XOR, YOR.XMX.MASK.fi, 0 ) 

C ALL  PLOTF ( 0 i i i ,XX , Y Y f ZZ , XSCALE , YSCALE > ZSCALE ,NN , NX , -PHI , 

*  -THETA. XOR, YOR.XMX. MASK. 0.0) 

41  CONTINUE 
100  RETURN 

END 

BLOCK  DATA 

COHMON/DATAi/XMARG . XWMAX . XINF , Y INF , XLINF . YWMAX . X YMAR . ASPR 
COMM0N/DATA2/IPLX.  IPLZ.ThETA. PHI .DELTA. XSC.YSC.ZSC. SIZE 
COMHON/DATA3/IPL P , INFORM .  ITITLE, IPERIM . IAXIS , IMASK 
C0MM0N/DATA4/XMIN , XMAX , ZMIN , ZMAX , YMIN . YM AX 
COMMQN/DATAS/NXPL.NZPl ,  IPL.PX , IPLPZ , XSCALO  .YSCALO  .ZSCALO 
COMMQN/DAT  A7/NHEDX . HEDX ( 5 ) . NHEDZ , HEDZ ( 5 ) ,NhEDY , HEDY ( S ) 
CnMM0N/DATA8/NTITLi jTITLEl (2) ,NTITL2.TITLE2(15) ,NTITL3.TITlE3( 15 > 


DATA  HEDX/4H 

,-4H 

,4H 

,4H 

>4H 

l 

/ 

DATA  HEDY/4H 

,4H 

.4H 

,4H 

,4H 

/ 

DATA  HEDZ/4H 

,4H, 

,4H 

.  4H 

.4H 

/ 

DATA  TITLE1/4HFIGU.4HSE  / 


DATA 

TITLES 

/4h 

.4H 

,4H 

,4H 

,4H 

.4H 

,4H 

C4H 

.  4H 

,4H 

,4H 

,4H 

,4h 

,4H 

,4H 

/ 

DATA 

TITL.E3 

/4H 

,4H 

,4H 

,4H 

,4H 

,4H 

,4H 

C4H 

,4H 

,4H 

,4h 

.4H 

,4H 

,4H 

,  4H 

/ 

END 

SUBROUTINE  PLTIF(XUMIN.XUMAX,YVMIN,YUMAX) 
COMMON/DATAi/XMARG, XWMAX .XINF  . YINF.Xl.INF  .  YWMAX . XYMAR . ASPR 
COMMON/DATA2/  IPLX .  IPLZ .  THETA  .PHI ,  DELTA. XSf. .  YSC .  ZSC ,  SIZF 
C0MM0N/DATA3/IPLP , INFORM , ITITLE , IPERIM , 3 AXIS , IMASK 
COMKQN/DATA4/XMIN . XMAX , ZMIN .ZMAX .YMIN , YMAX 
COMMON/DAT AS/NXPL . NZPL ,  IPLPX , IPLPZ , XSCALO . YSCALO , ZSCALO 
COMMON/DAT A7/NHEDX .  HEDX  ( S ) ,  NHEDZ . HEDZ  ( S ) ,  NHF.DY , HEDY ( S ) 
DATA  AN0/2HN0,  AYE/2HYE/ 

YSP=0.03*YUMAX 


V 
k  , 


•: 

DXi=. iOfiXVMAX 

:• 

DX2=. 125*XUMAX 

■ 

DX3=.0?*XVMAX 

1 

X0=XVMIN 

. 

• 

Y0=YVMIN 

-  ’ 

XOR=XINF 

/ 

YOR=YINF 

X  =  XOR 

/ 

ft 

Y  =  YOR 

i  . 

r 

v'. 

,\ 

-■* 

,v 

L 


CALL  L£TTR(X,Y,7HX  MIN  =,0.0,7) 
Y=Y-YSP 

CALL  LETTR (X,Y,7HX  MAX  =,0.0,7) 
Y=Y-YSP 

CALL  L.ETTR C X , Y , 7HZ  MIN  =,0.0,7) 
Y=Y-YSP 

CALL  LETTR (X,Y,7HZ  MAX  =,0.0,7) 
Y=Y-YSP 

CALL  LETTR CX, Y,7HY  MIN  =.0.0,7) 
Y=Y-YSP 

CALL  LETTR (X,Y,7HY  MAX  =,0.0,7) 

X=X+DXi 

Y-YQR 

CALL  NUMBR(X,Y,XMIN, 0.0,2) 
Y=Y-YSP 

CALL  NUMBR(X,Y,XMAX,fl .0  ,2) 
Y=Y-YSP 

CALL  NUMBS ( X , Y , ZMIN ,  0  0 , 2 ) 
Y=Y-YSP 

CALL  NUMBR(X,Y.ZMAX,0.0,2) 
Y=Y-YSP 

CALL  NUMBRLX.Y.YMIN, 0.0.2) 
Y=Y-YSP 

CALL  NUMBS  <  X , Y , YMAX ,0.0,2) 

Xi=  33HXVMAX 
Y)=YVMAX 

CALL  BOXIN<XO ,Y0,Xi,Yi ,0 .0 ,0 .0) 
XOR=Xl+XHARG 
X=XOR 
Y=YOR 

CALL  LETTR (X.Y.7HTHETA  =.0.0,7) 
Y='r-YSP 

CALL  IETTR(X,Y,7H  PHI  =,0.0,7) 
Y=Y— YSP 

CALL  LETTR (X,Y,7H  SIZE  =,0.0,7) 
Y=Y-YSP 

CALL  LrTTP(X,Y,7H  XSC  =.0  0,7) 
Y=Y-YSP 

CALL  LETTR(X,Y,7H  YSC  =,0.0,7) 
Y=f-fSP 

CALL  l£TTR(X,Y.7H  ZSC  =,0.0.7) 

x=x+dx2 

Y-YOR 

CALL  NUMBR(X,Y, THETA, 0.0,1) 
Y=Y-YSP 


■*  *  ■**-  ■  ■  •  -  —  •  ■-  , 


.... 
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CALL  NUMBR(X,Y, PHI, 0.0,1) 

Y=Y-YSP 

CALL  NUMBR(X,Y, SIZE, 0.0,1) 

Y=Y-YSP 
XSCO=XSC*SIZE 
YSCO=YSC*SIZE 
ZSCO=ZSC*SIZE 
CALL  NUHBR(X,Y.XSCO,0. 0,1) 

Y=Y-YSP 

CALL  NUMBR(X,Y,YSCO, 0.0,1) 

Y=Y-YSP 

CALL  NUMBR(X,Y,ZSCO, 0.0,1) 

X2=X1+.25*XVMAX 

Y1=YVMAX 

CALL  BOXIN(Xi,YO,X2,Yi .0 .0,0.0) 

X0R=X2+XMARG 
YOR=YQR-.  OKYVMAX 
Y3P=0 . 844YUHAX 

50  CALL  ZALPH(38,38H  DO  YOU  WANT  NEW  TITLES  FOR  THE  AXFS  ) 
CALL  ZKYBDfl ,4.NANSW. ANSW) 

IF (ANSW.EQ. ANO)  GO  TO  100 
IF(ANSW.NE.AYE)  GO  TO  50 

CALL  ZALPH(3S,3SH  TITLE  FOR  THE  Y-AXIS  (20  Chars.)  ) 
CALL  ZKYBD(1,20,NHEDY,HEDY) 

CALL  ZALPH(3S,3SH  TITLE  FOR  THE  X-AXIS  (20  Chars.)  ) 
CALL  ZK  YBD ( 1 , 20 , NHEDX , HEDX ) 

CALL  ZALPH(3S,35H  TITLE  FOR  THE  Z-AXIS  (20  Chars.)  ) 
CALL  ZKYBD(i,20,NHEDZ.H£DZ) 

10  FORMAT (A2) 

20  FORMAT ( 20 A4) 

100  X  =  XOR 
Y  =  YOR 

CALL  LETTR(X,Y,4HY  =  ,0.0,4) 

Y=Y-YSP 

CALL  LETTR (  X,Y,4HX  =  ,0.0,4  ) 

Y=Y-YSP 

CALL  LETTR <X,Y,4HZ  =  ,0 .0,4) 

Y=Y-i.5*YSP 

CALL  LETTR (X ,Y ,21H  PLT3T,  Uer.  3D:04:84  ,0.0,21) 

X=X+DX3 

Y=YOR 

CALL  LETTR ( X ,  Y , HEDY( 1 ) , 0 . 0 , NHEDY ) 

Y=Y-YSP 

CALL  LETTR (  X,Y,HEDX(i) ,0 . 0 ,NHEDX) 

Y=Y-YSP 

CALL  LETTR (  X,Y,H£DZ(1>,0.0,NHEDZ> 

X0=X2 

X3=XVMAX 

Yi=YVMAX 

CALL  B0XIN(X2,YO,X3,Yi,Q .0 ,0 .0) 

RETURN 

END 

SUBROUTINE  DATEP(X0,Y0) 

DIMENSION  IDATE(iS) 


CALL  FTIHE(IDATE) 

CALL  LETTR < XO , YO , IDATE ( 6 ) f 0 . 0 . 20 ) 

RETURN 

END 

SUBROUT INE  REFAX(XC , YC , ZC , XNAX , YHAX , ZMAX ) 
COMMON/DATAl/XNARG , XWHAX , XINF , YINF , XLINF , YWHAX , XYMAR . ASPR 
COMMON /DATA2/IPLX , IPlZ .THETA , PHI .DELTA ,XSC , YSC ,ZSC .SIZE 
DATA  RD/57.295779513082/ 

DATA  HTO/ .01/ .HTH0/0 .001/ 

THETAR=THETA/RD 

PHIR=PHI/RD 

HT=HTO 

UHT^NHTO 

SINPH  =  SIN(  PHIR  ) 

COSPH  =  COS(  PHIR  ) 

SINTH  =  SIN(  ThETAR  ) 

COSTH  *  COS(  THETAR  ) 

XLMX=XMAX' 

YLMXaYMAX 

ZLMX=ZMAX 

SX=XLMX*XSC/AB3tXSC) 

SY=YLMXIY3C/AFS ( YSC ) 

B7"ZLMXIZSC/ABS  t  ZSC ) 

XAXX  =  SX  *  COSPri 
YAXX  =  SX  *  SINPH  *  SINTH 
XAXY  =  0 

YAXY  =  SY  *  COSTH 

XAXZ  =  SZ  *  SINPH 

YAXZ  =  -SZ  *  COSPH  *  SINTH 

X  =  XC 

Y  =  YC 

CALL  ZMOUEt  X  ,  Y  ) 

X  =  X  +  XAXX 

Y  -  Y  +  YAXX 

CAL.  ZDRAUt  X  .  Y  ) 

X-X+HT 

Y=Y-2.0*HT 

CALL  LETTRtX.Y.lHX. 0.0,1) 

X  =  XC 

Y  =  YC 

CALL  ZMOVEt  X  ,  Y  ) 

X  =  X  +  XAXY 

Y  =  Y  t  YAXY 

CALL  ZDRAWt  X  ,  Y  ) 

X-X+HT 

Y=Y-HT 

CALL  LETTR (X.Y.iHY. 0.0,1) 

X  =  XC 

Y  =  YC 

CALL  ZMQVEt  X  .  Y  ) 

X  =  X  ♦  XAXZ 

Y  a  y  +  YAXZ 

CALL  ZDRAWt  X  ,  Y  ) 


Y=Y-2 . *HT 

CALL  LETTR(X,Y,1HZ. 0.0,1) 


CALL  ZMOVEf  X  ,  Y  ) 

RETURN 

END 

SUBROUTINE  BOXIN ( XO .. YO , XMAX , YHAX , XYMAR , ASP R ) 

X«AR=ASPR*XYNAR 
YHAR-XYHAS 
Xi=X0+XMAR 
Yl=YO+YHAR 
X2=XHAX-XMAR 
Y?=YMAX-YMAR 
CALL  ZMOVE(Xi,Yi) 

CALL  ZDRAW(X1,Y2> 

CALL  ZDRAU(X2,Y2) 

CALL  ZDRAW(X2,Yi  ) 

CALL  ZDRAU(X1.Y1  ) 

RFTURN 

END 

SUBROUTINE  TITLE(XTMIN,XT«AX,YT«IN,YTMAX) 

COHhON/DATAi/XilASG ,  XvlMAX  ,  XINF ,  Y  INF ,  XLINF ,  YHHAX ,  XYMAR .  ASPR 
COMMON/DAT A2/IPLX ,IPLZ . THETA ,  PHI , DELTA , XSC , YSC , ZSC ,SI ZE 
C0HM0N/DATA4/XNIN .  XNAX .  ZHIN .  ZHAX ,  YillN ,  YMAX 
C0MN0N/DATA8/NTITL1 ,TITLE1(2) ,NTITL2.TITLE2(iS) ,NTITL3,TITLE3(i5) 
DATA  ANO./2HNO/AYE/2HYE/ 

YSP=0  S*YTMAX 

10  CALL  ZALPH(31 ,3iH  DO  YOU  WANT  NEW  FIGURE  TITLE  ) 

CALL  ZK  YBD 1 1 , 4 , NANSW . ANSU ) 

IF (ANSW.EQ.ANO)  GO  TO  100 
IF'. ANSU. NE. AYE)  GO  TO  10 

CALL  ZALPH(3i .31HTYPE  THE  1ST  LINE  (60  Chars.)  > 

CALL  ZXYBD(1 .60  ,NTITL2,TITLE2) 

CALL  ZALPH(3i .31HTYPE  THE  END  LINE  (60  Chars.)  .) 

CALL  ZK YBD( 1 .60 .NTITL3 .TITLE3) 

100  X0=XTHIN 
Y0=YTMIN 
XQR=X0+XMARG 
YOR=YTKAX-Xf!ARG 
X=XOR 
Y=YOR 

CALL  LETTR(X,Y,TITLE2(1),0  0.NTITL2) 

Y=Y-YSF 

CALL  LETTR (X , Y ,TITLE3( 1 ) , 0 . 0 .NTITL3 ) 

Y=Y-YSP 

C ALL  BOX I N ( X0 , Y  0 , XTMAX . YTMAX . 0 . 0 , 0 . 0 ) 

RETURN 

END 

SUBROUTINE  DEVON 
DINENSION  SUBR ( 3 ) 

DATA  SUBR/4H  .4H  ,4H  / 

LU=i 


CALL  7BEGN 


CALL  ZAINTd  ,  IERROR) 

IFdERROR  .EQ.O)  GO  TO  100 
CALL  ERRORS! 1 , 1ERROR , 6HZAINT  ) 

C  GO  TO  999 
100  CALL  ZK INT  < 1 , TERROR ) 

IFdERROR. EQ.O)  GO  TO  200 
CALL  ERRORS ( 1 , IERROR , 6HZK INT  ) 

C  GO  TO  999 
200  CALL  ZBINTd,  IERROR) 

IF (IERROR. EQ.O)  GO  TO  400 
CALL  ERRORSd, IERROR, 6HZBINT  ) 

C  GO  TO  999 
400  ICOTR=0 

CALL  ZDINT(LU,ICOTR, IERROR) 

IF < IERROR. EQ  0)  GO  TO  600 
CALL  ERRORSd,  IERROR, 6HZDINT  ) 

C  GO  TO  999 
600  CONTINUE 
999  RETURN 
END 

SUBROUTINE  ERRORS(LU , IERROR ,SUBR ) 

DIMENSION  SUBR(3) 

CALL  ZMCUR 

URITE(LU,100)  IERROR, SUBR 

100  FORHAT("<<<<<  Error  ".IS,"  occured  in  subroutine”. 3A4, “>>>>>") 
RETURN 
END 

SUBROUTINE  VIEUP (XUMAX , YUMAX ,SIZE , IERROR ) 

DIMENSION  AR(2) ,UIEU(4) 

CALL  ZIWS(2S4,0, 2, IDUM.AR, IERROR) 

XSIZE=i.fl 

YSIZE=AR(2) 

CALL  ZASPK(XSIZE,YSIZE> 

WXMIN=0 . 0 
UXMAX=XUMAX 
UYMIN=0 . 00 
UYMAX=YUMAX 

CALL  ZUIND(tiXMIN,UXMAX,UYMIN,WYMAX) 

CALL  ZIUS ( 451 , 0 , 4 , IDUM . VIEU , IERROR ) 

YMINX=.0S*VIEU(2)/SIZE 
VHAXX*.9SWIEH<2)/8IZE 
VMINY=.0S*UIEU(4)/SIZE 
VMAXY=.95*VIEU(4)/SIZE 
CALL.  ZVIEN(VMINX,WIAXX,VMINY,VMAXY) 

XCSIZ=0 . 0 15# ( WXMAX— UXMIN ) /SIZE 
YCSIZ=0 . 02S*(WYMAX-WYMIN)/SIZE 
CALL  ZCSIZ(XCSIZ,YCSIZ) 

RFTURN 
END 

SUBROUTINE  VPMAX(UIEW) 

DIMENSION  AR(2),9IEM(4) 

CALL  ZIWS(4Si, 0,4, IDUM, VIEW, IERROR) 

CALL  ZIQS(254 ,0,2, IDUM, AR , IERROR ) 

XNEM-i.O 
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YNEW=i.O 

IF(AR(2)  LE . 1 .0)  YNEW=AR(2) 

IF ( AR < 2 )  GT.i.O)  XNEU=i./AR<2) 

CALL  ZVIEW(0 . 0 , XNEU ,0.0 , YNEW) 

RETURN 

END 

SUBROUTINE  DEVQF 
CALL  ZCOLR(O) 

CALL  ZAEND 
CALL  ZBEND 
CALL  ZKEND 
CALL  ZDEND 
CALL  ZLEND 
CALL  ZEND 
RFTURN 
END 

SUBROUTINE  LETTR  (X,Y,IAR,ANG,NC) 

DIMENSION  IAR(i) .VIEw(4) 

CALL  ZMOVE(X,Y) 

CALL  VPMAX(VIEW) 

CALL  ZTEXT(NC.IAR) 

CALL  ZVIEW (VIEW C i ) ,UIEW(2) ,UIEN(3) ,V1EW<4) ) 

RETURN 

END 

SUBROUTINE  NUMBR ( X , Y , FPN , THETA , IFMT ) 
DIMENSION  JFMT (S) ,ISYM(3i 
GO  TO  (10,20,30)  IFMT 
10  ENCODE  US,  IS.  ISYH)  FPN 
GO  TO  100 

20  ENC0DE(15,2S,ISYM)  FPN 
GO  TO  100 

30  ENCODE* 10 ,35 ,ISYM>  FPN 
LOO  CALL  LF.TTR(X,Y.ISYM, THETA. 15) 

15  FORMAT (F5.1) 


FORMAT (F5.0) 

RETURN 

END 
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FTN7X  .0 
OFILES  0,1 
$ENA/BLKi/ 

PROGRAM  PLT3P (S) .  G.A.KERAMIDAS  3P  =  04  =  25  =  84 
COMMON/BLKi/MASK(9000)  ,XPLd51) .ZPL(iSi) .YPLdSi  .151) 
COMMON/OATA1/XMARG , XWMAX , XINF . YINF, XL INF , YWMAX , XYMAR .  ASPR 
COMMON/DAT A2/IPLX, IPLZ, THETA , PHI  .DELTA , XSC ,  YSC , ZSC  ,SIZF. 
COMMON/DAT A3/IPLP , INFORM . ITITLE ,  IPERIM , IAXIS , IMASX 
C0NM0N/DATA4/XMIN . XMAX ,ZMIN, ZMAX , YMIN , YNAX 
COMHON/DATA5/NXPL -NZPL , IPLPX , IPLPZ , XSCALC . YSCALD , ZSCALO 
COMMON/DAT A6/INST  N , IFAST , IPLOT , IPMENU . IPEN< 2) 

C0MM0N/DATA7/NHEDX . HEDX ( 5 ) , NHEDZ , HEDZ < S  > , NHEDY , HEDY ( S ) 

COMMON/DAT A8/NTITL1 . TITLE! ( 2 ) .NTITL2. TITLES (lb) .NTITL3.TITLE3(15) 
DIMENSION  VIEW ( 4 ) 

Y0=0  0 
X0=0.0 
CALL  DEVON 

CALL  VIEW  (XWMAX ,  YWMAX ,  SIZE ,  I  ERROR ) 

CALL  ZIUS (451 , 0 , 4 . IDUM , VIEW , IERROR ) 

IFdERROR  EO.fl)  GO  TO  305 

CALL  ERRU5U.  IERROR,  1CHZIWS.  45 1  > 

GO  TO  999 

305  IF ( INFORM. LE.O)  GO  TO  310 
Yi=0.8*YWMAX 
YPMAX=.2*VIEW<4) 

CALL  PLTIF<X0, XWMAX, Yi, YWMAX) 

310  IF< ITITLE. LE.O)  GO  TO  315 
Y2=0.08*YWMAX 

CALL  TITLE ( XO , XWMAX, Yfl ,Y2  > 

YQ=.08*VIEW<4> 

315  CONTINUE 

CALL  ZIWS(45i ,0,4, IDUM .VIEW , IERROR ) 

IFdERROR. E9.0)  GO  TO  400 

CALL  FRRORSd, IERROR,  10HZ1US.  451  ) 

GO  TO  999 
400  XV0=VIEW(i) 

YV0=VIEW(3)+Y0 

XVMAX=VIEW(2) 

YVMAX=VIEW(4)-YPMAX 
410  CALL  ZVIEU(XVO.XVMAX,YVO,YVMAX) 

CALL  ADISP 

CALL  ZALPH  <40,40H  ****#*PLOTTING . !  !!!«*«*«  ) 

CALL  ZCOLR(l) 

CALL  BOXIN ( . 0 , . 0 , XWMAX , YWMAX ,00,0.0) 

CALL  ZC0LRdPEN<2>> 

C ALL  P LOTS ( M ASK , XPL , ZPL , YP L , NXP  L , NZP  L .  1 ,  NXP L ) 

CALL  ZMCUR 
500  CONTINUE 
999  CONTINUE 
CALL  DEVQF 
CALL  SEGRT  (Z991) 

END 

SUBROUTINE  ADISP 

INTEGER  ILIST <7) ,  STRINGS),  IERR 
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DATA  STRING  /15550B.1SS12B/ 


CALL  ZALPH  (4. STRING) 

RETURN 

END 

SUBROUTINE  PLOTS  <  MASK .  X ,  7. ,  Y ,  NX  ,NZ ,  IDI H ,  I  DIHEN ) 

LOGICAL  LPLOTX, LPLOTZ 

DIMENSION  XX(i5i),ZZ(i5i),YY(i5i),U£RTEX(1.6) 

EHA  MASK (9BQQ).X( 151) (Z(i5i) ,Y(lSl ,151) 
COMMQN/DATAi/XHARG.XNMAX,XINF,YINf ,XLINF , YWMAX , XYKAR ,ASPR 
C0MM0N/DATA2/IPLX .  IPLZ , THETA , PHI , DELTA . XSC , YSC , ZSC  SI 7E 
COMM0N/DATA3/IPLP .INFORM , ITITLE , IPERIM , I AXIS , IMASK 
C0MN0N/DATA5/NXPl . NZPL , IPLPX. IPLPZ . XSCALO , YSCALO , ZSCAL  0 
COMMON/DATA4/XMIN , XHAX , ZMIN . ZMAX , YHI N . YHAX 
DATA  SU/  -1 . / 

PID180=. 01745329252 
PIR=3. 14159 

COSPH  =  COS (  PID180  *  PHI  > 

COSTH  =  COS(  PID180  *  THETA  i 
STNPH  =  SIN(  PID130  *  PHI  ) 

SINTK  =  SIN (  PID180  *  THETA  ) 

TMP1  =  XSC  *  COSPH 

TMP2  =  ZSC  *  SINPH 

XMX  =  AMAXK  TMP1  ,  TMP2  ,  TM?1+TMP2  ,  0.  ) 

XHN  =  AMINK  THP1  ,  TMP2  ,  TMP1+TMP2  ,  0.  ) 

TMP1  =  XSC  *  SINPH  KSINTH 

TMP2  -  -ZSC  4C0SPH  4SINTH 

YMX  =  AMAXK  TMPi  .  TM?2  ,  TMP1+TMP2  .  0.  ) 

YMN  =  AHINK  TMPi  ,  TMP2  ,  TMP1+TMP2  .  0.  > 

TMPI  =  YSC  *  COSTH 

YMX  =  AMAXK  YMX  ,  YMX  +  TMPI  > 

YMN  =  AMINK  YMN  ,  YMN  +  TMPi  > 

XL  =  XMX  -  XMN 
Yl.  =  YMX  -  YMN 

IF<  XL  .GT.  XWMAX  OR.  YL  .GT.  YNMAX  )  CO  TO  998 

XOR  =  XMARG  -  XMN 

XMX  XMARG  +  XL  +  XMARG 

S'-I=-SU 

YOR  =  ABS ( YMN ) + ( YHM AX-YL ) /2 . +DELT A 
LPLOTX  =  -SINPH  *  COSTH  GT.  0. 

LPLOTZ  =  COSPH  *  COSTH  .GT.  0. 

XSCALE  =  XSC/(XMAX-XMIN> 

Y SCALE  =  YSC/(YMAX-YMIN) 

ZSCALE  =  ZSC/(ZMAX-ZMIN) 

CALL  ZMQVE(  XCR  .  0.0) 

KPLX  =  MAX0(  IPLX  ,  1  ) 

NEXTX  =  1 

DO  450  NLINE  -  1 ,NX 

IF(  NLINE  .NE.  NEXTX  )  GO  TO  450 

NEXTX  =  MIN0(  NLINE  +  KPLX  .  NX  ) 

NN  =  NLINE 

IF(  LPLOTX  )  NN  =  NX  +  1  -  NLINE 
DO  400  J=i  ,  NZ 


XX(J)  =  X (NN)-XMIN 
YY(J)  =  Y (NN, J)-YMIN 
Z7(J>  =  Z(J)-ZMIN 
400  CONTINUE 
NN  =  NLINE 

CALL  P LOTF (  i i i i , XX , YY , ZZ , XSCALE , YSC ALE . ZSCALE , NN , NZ , -PHI , -THETA 

*  XOR, YOR, XHX, MASK, VERTEX) 

4S0  CONTINUE 

IF(IPERIM  GE.i) 

*C ALL  PER IM<  X , Y , Z ,XX , YY , ZZ , XSCALE , YSC ALE . ZSCALE ,XHIN , YHIN . ZMIN , 
iNX , NZ . PHI , THETA , XOR , YOR , XMX , MASK . NN , i ) 

KPLZ  =  MAX0<  IPLZ  ,  i  ) 

NEXTZ  =  i 

DO  SSO  NLINE  =  i  ,  NZ 

IF(  NLINE  .NE.  NEXTZ  )  GO  TO  SSO 

NEXTZ  =  HINO (  NLINE  +  KPLZ  ,  NZ  ) 

NN  =  NLINE 

IF(  LPLOTZ  )  NN  =  NZ  +  i  -  NlINE 
DO  500  I  =  i  .  NX 
XX(I)  =  XCD-XMIN 
YY(I)  =  Y(I,NN)-YNIN 
mi)  -  Z (NN)-ZMIN 
500  CONTINUE 
NN  =  NLINE 

CALL  P LOTF (  0 i i i , XX , YY , ZZ , XSCALE , YSC ALE , ZSCALE , NN , NX . -PHI , -THETA 

*  XOR , YOR , XHX , MASK .00) 

SSO  CONTINUE 

IF(IPERIM.GE.i) 

ICALL  PERIH(X,Y,Z, XX, YY,ZZ, XSCALE, YSCALE, ZSCALE, XMIN.YMIN, ZMIN, 
INX , NZ , PHI ,  THETA , XOR , YOR , XHX , MASK , NN , 2 ) 

XC0=XWMAX-4.*XMARG 
YC0=  8SIYWMAX 
ZC0=XCC 

XLMX=2.S*XMARG 
YLMX=3.fl#XHARG 
ZLMX=2. 54XMARG 
CALL  ZCOLR(l) 

IF(IAXIS.GE.l)  CALL  REFAX(XCO,YCO,ZCO,XLHX,YLMX,ZLMX> 

XP0=XMARG 
YPG=.95*YWMAX 
CALL  DATEP (XPC . YPO ) 

GO  TO  999 
998  CONTINUE 

CALL  LETTR (0.2,0 .25.  02,2bHFIGURE  SIZE  LIMIT  EXCEEDED, 0 . ,26  ) 

999  CONTINUE 
RETURN 
END 

SUPRGUTINE  PLOTF< IVXYZ , XDATA , YDATA ,ZDATA , XSCALE , 

JtYSCALE ,  ZSCALE  .NLINE ,  NPNTS ,  PHI .  THETA ,  XREF , 

IYPEF, XLENTH , MASK , VERTEX ) 

INTEGER  HIGH,  QLDHI ,  CLDLOW 

COMMON/DATA2/IPLXO,:PLZO,THETAO,PHIQ, DELTAS, XSCO,YSCO,ZSCO, SIZE 
COMMON/DATA3/IPLP , INFORM, ITITLE , IPERIM, IAXI3, IMASK 
COMMON/DAT  A5/NXPL , NZPL , IPLPX . IPLPZ , XSCALO , YSCALO , ZSCALO 


DIMENSION  XDATA(i),  YDATA(l),  ZDATA(i) ,VERTEX(16) 

EMA  MASK (9000) 

DATA  INIT,  JUXYZ,  SPHI,  STHETA/-i,-i,-1.0E30,-1.0E30/ 
IF  (NLINE.EQ.O)  GO  TO  550 
IF  (NLINE.NE.i)  GO  TO  20 
PIPI=iOO . (SQRT (SIZE) 

NYPI=PIPI*20 

LIMITX= ( XLENTHIP IP  I ) +0 . 5 
I=LIMITX+LIMITX 
DO  10  K=i,I 
MASK(K)=INIT 
I  CONTINUE 
INIT=-i 
INCI=-i 
1=0 

20  IF(JVXYZ.EQ.IVXYZ)GO  TO  70 
JVXYZ=IUXYZ 
INDZ=i 
INDY=i 
INDX=i 
INDV=1 

X F ( JV) X YZ .  LT  .1000) GO  TO  30 
INDV=2 

JVXYZ=JVXYZ-1000 
30  IF(JVXY2.LT  iOO)GQ  TO  40 
INDX=2 

JVXYZ=JVXYZ-100 
40  IF(JVXYZ.LT. 10)GO  TO  50 
INDY=2 

JVXYZ=JUXYZ-10 
50  IF  ( JUXYZ  LT.DGG  TO  60 
INDZ=2 

60  JUXYZ=IVXYZ 

70  IF  ( PHI. EQ. SPHI. AND. THETA. EO.STHETA)  GO  TO  80 
SPHI=SIN(.017453292S*PhI) 

CPHI=COS(.0174532925*PHI) 

STHETA=SIN(  01745329254THETA) 

CTHETA=COS( . 0I7453292SITHETA) 

Aii=CPHI 
A)3=-SPHI 
A21=STHETA*SPHI 
A22=CTHETA 
A23=STHETA*CPHI 
SPHI=PHI 
STHETA=THETA 
80  INCI—INCI 

IF(I.NE.0)I=NPNTS+1 
DO  530  K»i ,NPNTS 
I=I+INCI 

GO  TO  (90,100), INDX 
90  X=XDATA( 1 ) + ( 1-1 ) XXSCALE 
GO  TO  110 

100  X=XDATA(  I XXSCALE 
110  GO  TO  (120, 130), INDY 


120  Y=YDATA(1 )+(I-l )*YSCALE 
GO  TO  140 

130  Y=YDATA(I)*YSCALE 
140  GO  TO  (150,160) ,INDZ 
150  Z=ZDATA<l)+(NLINE-i)*ZSCALE 
GO  TO  170 

160  Z=ZDATA(I)*ZSCALE 
170  XXX=A11*X+A13*Z+XREF 
XX=XXX 

IX=(XX*PIPI)+.5 
YYY=A2i*X  +  A23IZ  +  YREF 
YY=YYY+A22*Y 
IY=(YY*PIPI)+.5 
IF  (IX.LE.0)IX=1 
IF ( IX . GT . LIMITX) IX=LIHITX 
IF(IY.LT.5)IY=5 
IF ( I Y . GT . NYPI ) I Y=NYP I 
IF  (K.NE.l)GO  TO  250 
LH'J=IX+IX 
HIGH  =LOW-i 
MLOU=MA5X(LOW> 

MHIGH=MASX (HIGH) 

IF(MHlGH-IY)2Gfl,210 ,180 
180  IF(MLOU'-IY)  190 .230,220 
190  LOCQLD=G 
GOTO  240 

200  MASK ( HIGH ) =1 Y 

IF(MLOW.EQ.-l) MASK ( LOW ) = I Y 
210  lOCOLD=+1 
GO  TO  240 

220  MASK  <  LOW i=I Y 

231)  LQCQLD=-1 

240  IF ( IMASK )  241,241,242 

241  CALL  ZMOOE(XX,YY) 

GO  TO  243 

242  IF ( LOCOLD . NE . 0 )  CALL  ZMCVE  < FLOAT ( IX ) /P IP  I , FLOAT ( I Y ) /P IP I ) 

243  CONTINUE 
JX=IX 
JY=IY 
IYREF=IY 

IFdNDV.EQ.l )GO  TO  530 
INDEX=INCI+6 
VERTEX! INDEX)=XX 
VERTEX ( INDEX+i )=Y Y 
VERTEX  < INDEX+8 )=XXX 
VERTEX ( INDEX+9 )= Y  YY 
IF (NLINE.NE.l)  GO  TO  530 
VERTEX ( 1 )=XX 
VERTEX ( 2 )=YY 
VERTEX  <  9 ) =XXX 
VERTEX(iO)=YYY 
GO  TO  530 

250  IFdMASK .GT. 0)  GO  TO  255 
CALL  ZDRAW(XX.YY) 
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GO  TO  530 

255  IFdX.NE.  JX)GO  TO  260 
JY=IY 
GO  TO  280 

260  Y INC=FLOAT ( I Y- J Y ) / ABS ( FLOAT ( I X- JX ) > 

INCX= ( IX-JX ) /I ABS ( IX-JX ) 

YJ=JY 

270  JX=JX+IHCX 
YJ=YJ+YINC 
JY=YJ+.S 
L0U=JX+JX 
HIGH=L0U-i 
NL0U=MASK(LOW) 

MHIGH=MASK(HIGH) 

280  IF(MHIGH-JY)3fl0 ,300 ,290 
290  IF(HLOW-JY>310,320,320 
300  LOC=+i 

IF (LOCQLD)360 ,370 ,430 
310  LQC=0 

IF( LOCOLD) 340 , 350 , 330 
320  LOC=-i 

I F ( LOCOLD  >510 ,450,440 

330  IF (MHIGH . LE . IYREF ) CALL  ZDR AU < FLOAT  <  JX ) /P IP  I , FLOAT  (MHIGH > /P IP  I > 
GO  TO  350 

340  IF(MLOU .GE . IYREF )CALL  ZDR AU(FL0AT ( JX) /P IP I , FLOAT ( MLOU ) /PIP I ) 
350  CONTINUE 
GO  TO  520 

360  IF (HLQU-IYREF)370 ,380,380 

370  IF(HHIGH-IYREF)400 ,390,390 

380  CALL  ZDR AU ( FLOAT  <  JX ) /P IP I , FLOAT ( MLOU ) /P I P I ) 

390  CALL  ZMOVE (FLOAT ( JX ) /P IP I , FLOAT (MHIGH ) /P IP  I ) 

GO  TO  430 

400  IFtHHIGH.EQ  -l)GO  TO  430 
OLDHI =HI GH-20I NCX 
IF < MASK < OLDHI )-JY> 420, 420, 410 
410  CALL  ZMOVE ( FLOAT  <  JX ) /P IP  I , FLOAT  <  JY ) /P IP  I ) 

GO  TO  430 
420  IP=JX-INCX 

ID=NASK (OLDHI) 

CALL  ZHOVE ( FLOAT ( IP ) /P IP I , FLOAT ( ID ) /P IP I ) 

430  MASK < HIGH) =JY 

IF(ML0U . EP . -1 )MASK (LQU)=JY 
CALL  ZDRAU( FLOAT ( JX ) /P IP I , FLOAT ( J Y ) /P IP I ) 

GO  TO  520 

440  IF (MHIGH-IYREF >460,460, 450 

450  IF (MLOU-IYREF ) 470 ,470,480 

460  CALL  ZDR AU ( FLOAT <  JX ) /P IP I , FLOAT  <  HH IGH ) /P IP I ) 

470  CALL  ZMOVE ( FLOAT ( JX ) /P I P I , FLOAT ( MLOU  > /P I P I ) 

GO  TO  510 

480  0LDL0U-L0U-2*INCX 

IF(MASK(OLDLOW)-JY)490, 500,500 
490  CALL  ZMOVE (FLOAT (JX)/PIPI,FLOAT(JY)/PIPI) 

GO  TO  510 
500  IPsJX-INCX 
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I  D=H  ASK  <  OLDLOuJ ) 

CALL  ZMQVE < FLOAT <  IP  >/P  IP  I , FLOAT (I D > /P IP  I ) 

510  MASK ( LOW )=JY 

CALL  ZDRAW ( FLOAT ( JX ) /? IP I , FLOAT ( JY ) /P IP I ) 

520  IYREF=JY 
LOCOLD=LOC 

IF(JX.NE.IX)  GO  TO  270 

530  CONTINUE 

IF(IHASK)  531,531,532 

531  CALL  ZMOVE(XX,YY) 

GO  TO  535 

532  CALL  ZMOVE ( FLOAT ( JX ) /P IP  I , FLOAT  <  JY ) /P IP  I > 

535  IFdNDV.EQ.DGO  TO  540 

lNDEX=-INCI+6 
VERTEX (INDEX )=XX 
VERTEX ( INDEX+i )=YY 
VERTEX* INDEX+8)=XXX 
VERTEX (INDEX+?)=YYY 
IF  (NLINE.NE.DGO  TO  540 
VERTEX ( 3 )=XX 
VERTEX(4)=YY 
V£RTEX(11)=XXX 
VERTEX! 12)=YYY 
540  1=1-1 
RETURN 
55 G  INIT=0 
RETURN 
END 

SUBROUTINE  PERIM^Y^XX.YYjZZjXSCALEjYSCALEjZSCALEjXMI^YMIN, 
*  ZMIN, NX , NZ , PHI , THETA , XOR , YOR , XMX , MASK .NN . IDIM ) 

COMMON /DATA3/IPLP , INFORM, ITITLE , IPERIM , IAXIS , IMASK 
C0MM0N/DATA5/NXPL , NZPL , IPLPX , IPLPZ , XSCAL0 , YSCAL0 , ZSCALO 
DIMENSION  XX(1) ,YY(1) ,ZZ(1) 

C0MM0N/DATA6/INSTN , IFAST , IPLOT , IPMENU , IPEN( 2) 

DIMENSION  XV(2) ,YV(2) ,ZV(2) 

EMA  MASK (90 00)  ..X(15i)  ,Z<i51),Y<15i,15i) 

IDIREC=-1 

IMIPLPX.GT.O)  GO  TO  2 
IPLPX=NX-i 
IPLPZ=NZ-i 
2  GO  TO  <4,ii),IDTM 

4  DO  5  J-l ,NZ 
XX(J)=X(1)-XMIN 
YY(J)=0.fl 

ZZ( J)=Z(J )-ZhIN 

5  CONTINUE 

DO  10  J=1,NZ, IPLPZ 

IDIREC=-ID1REC 

XV(i)=XX(J) 

XV(2)=XV(i> 

YV(l)=0.fl 
YV(2)-0 . 0 

IF(IDIREC.EO.-l)  YV(i)=Y(l,J)-YHIN 
IFdDIREC.Efl.+i )  YV(2)=Y(i , J)-YMTN 
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ZV(i)=2Z( J) 

ZV(2)=ZV(i) 

10  CALL  PLOTF<Oiii, XV, YU, ZO, XSCALE, YSCALE, ZSCALE, NN, 2, -PHI, 

*  -THETA, XOR, YOR, XMX, HASK, 0.0) 

CALL  P LOTF < fltii , XX , YY , ZZ , XSCALE , Y SC ALE , ZSCALE ,NN,NZ,-PHI, 

*  -THETA , XOR , YOR , XMX , MASK ,0.0) 

GO  TO  21 

11  DO  IS  1=1, NX 
XX(I)=X(I)-XMIN 
YY<I)=0.0 
ZZ(1)=Z(NZ)-ZMIN 

IS  CONTINUE 

DO  20  1=1 ,NX,IPLPX 

IDIREC=-IDIREC 

XV(i)=XX(I) 

Xy(2)=XV(i) 

YV(1)=0  0 
YV(2>=0.3 

IF(IDIREC.EQ.-l)  YV (i)=Y(I,N7)-YMIN 
IF(IDIREC.EQ  +1)  YU(2)=Y(I,NZ)-YMIN 
ZV<i)=ZZ(I) 

ZV(2)=ZV(i) 

20  CALL  P LOTF (0111, XV , Y V , ZW , XSCALE , YSCALE , ZSCALE , NN , 2 , -PHI , 

*  -THET A , XOR , YOR , XMX , MASK ,0.0) 

C ALL  P LOTF ( 0 1 1 1 , XX , YY , ZZ , XSCALE , YSCALE , ZSCALE , NN , NX , -PHI , 

*  -THETA , XOR ,  YOR , XMX , MASK ,0.0) 

GO  TO  31 

21  DO  25  J=1,NZ 
JJ=NZ+i-J 
XX(J)=X<NX)-XMIN 
YY(J)=0.0 
ZZ(J)=Z(JJ)-ZNIN 

2S  CONTINUE 

DO  30  J=1,NZ,1PLPZ 
IDIREC=-IDIREC 
JJ=NZ+1-J 
XV(i)=XX(J) 

XV(2)=XV(1) 

YV<l)-=0.3 
YV(2)=0 . 0 

IF(IDIREC.EO.-l)  YV(1)=Y<NX,JJ)-YMIN 
IF(IDIREC.EQ.+1)  YV(2)=Y(NX,JJ>-YMIN 
7V(1)=ZZ<J) 

zu(2)=zy<i) 

30  CALL  PLOTF ( 0 1 1 1 , XU , YU , ZV , XSCALE , YSCALE , ZSCALE , NN , 2 , -PHI , 

*  -THETA , XOR . YOR , XMX , MASK ,0.0) 

CALL  PLOTF ( E 1 11 , XX , YY ,ZZ , XSCALE , YSCALE , ZSCALE ,NN,NZ , -PHI , 

*  -THETA . XDP , YOR , XMX , MASK ,0.0) 

CO  TC  41 

31  DO  3S  1=1. NX 
II-NX+l-I 
xx(D=X(ii:-x«iH 
YY(I)=0.0 
ZZ(I)=Z(l)-ZMIN 

40 
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35  CONTINUE 

DO  40  I=i,NX,IPLPX 
II=NX+1-I 
IDIREC=-IDIREC 
XV(i)=XX(I) 

XV<2)=XV(i) 

Yg<i)=o  o 
YU(2)=YV(i) 

IF(IDIREC.EQ.-i)  YU(1)=Y<II,1 )-YMIN 
IF ( IDIREC . EQ . +i )  YV<2 )=Y ( I I , i >-YMIN 
ZV(1)=ZZ(I) 

zu<2>=zu(i> 

40  CALL  PLOTF <  0 1 1 1 ,  XU , YV , ZV , XSC ALE , Y SCALE , ZSC ALE ,  NN , 2 . -PHI , 

*  -THETA , XCR , YOR , XMX , MASK ,0.0) 

CALL  PLOTF ( 0 i i i , XX , Y Y , ZZ , XSCALE , Y SCALE , ZSCALE , NN , NX . -P HI , 

*  -THETA, XQR, YOR, XMX, MASK, 0 . G) 

41  CONTINUE 
100  RETURN 

END 

SUBROUTINE  PLTIF ( XUMIN . XUMAX , YUMIN , YUMAX ) 

COMMON/DATA1 /XMARG , XWHAX , XINF , YINF . XLINF , YWMAX , XYMAR . ASPR 
COMMON/DATA2/IPLX , IPLZ , THETA , PHI , DELTA , XSC , YSC , ZSC , SIZE 
C0MM0N/DATA3/IPLP , INFORM, ITITLE . IPERIrt , IAXIS , If ASK ' 
COMMON/DATA4/XMIN , XMAX , ZMIN , ZHAX , YMIN , YMAX 
COMMON/DAT A5/ NXPL , NZPL , IPLPX , IPLPZ , XSCALO , YSCALO , ZSCALO 
C0MM0N/DATA6/INSTN , IFAST , IPLOT , IPMENU . IPEN<2) 
COMMON/DATA7/NHEDX , HEDX  <  5 ) , NHEDZ , HEDZ ( 5 ) , NHEDY , HEDY ( 5 ) 
DATA  AN0/2HN0/AYE/2HYE/ 

YSP=0 . 034YVMAX 
DX1=.10*XWAX 
DX2=  12S*XUMAX 
D>3=. 074XUMAX 
X0=XVMIN 
YO-Y'.'MIN 
XOR=X‘NF 
Yoa»riNF 
X  =  XOR 
Y  =  YOR 

CALL  LETTR(X,Y,7HX  MIN  =,0.0,7) 

Y=Y— YSP 

CALL  LETTR(X,Y,7HX  MAX  =,0.u,7) 

Y=Y-YSP 

CALL  LETTR (X , Y ,7HZ  MIN  =,0  0,7) 

Y=Y-YSP 

CALL  LETTR(X,Y,7HZ  MAX  =,0.0,7) 

Y=Y-YSP 

CALL  LETTR (X , Y ,7HY  MIN  =  ,C.0,7) 

Y=Y-YSP 

CALL  LETTRIX,Y,7HY  MAX  =,0.0,7) 

X=X+DXi 


Y=YOR 

CALL  NUMBR(X,Y,XMIN,0 .0,2) 
Y=Y-YSP 

CALL  NUMBR<X,Y, XMAX, 0.0,2) 


Y=Y-YSP 

CALL  NUHBR(X,Y,ZMIN,0 .0,2) 
Y=Y-YSP 

CALL  NUHBR(X,Y,ZHAX,0 .0,2) 
Y=Y-YSP 

CALL  NUHBWX, Y,YrtIN,0. 0,2) 
Y=Y-YSP 

CALL  NUHBR(X,Y,YMAX,0 .0,2) 
X1=.33*XVMAX 
Yi=YVMAX 

CALL  BQXIN(XO,Yfl,Xl,Yl,0 .0,0. 0) 
X0R=X1+XHARG 
X=XOR 
Y=YOR 

CALL  LETTR(X,Y,7HTHETA  =,0.0,7) 
Y=Y-YSP 

CALL  LETTR(X,Y,7H  PHI  =.0.0,7) 
Y=Y-YSP 

CALL  LETTR(X,Y,7H  SIZE  =,0.0,7) 
Y=Y-YSP 

CALL  LETTR(X,Y,7H  XSC  =,0.0,7) 
Y=Y-YSP 

CALL.  LETTR(X,Y,7H  YSC  =,0.0,7) 
Y=Y-YSP 

CALL  LETTR(X,Y,7H  ZSC  =,0.0,7) 
X=X+DX2 
Y=YOR 

CALL  NUH8R(X,Y, THETA, C.0,1) 
Y=Y-YSP 

CALL  NUMBR(X,Y, PHI, 0.0,1) 
Y=Y-YSP 

CALL  NUhBRLX.Y, SIZE, 0.0,1) 
Y=Y-YSP 
XSCQ=XSC*SIZE 
YSCO=YSC*SIZt 
7SC0=ZSC*SIZE 
CALL  NUMBR(X,Y,XSCO, 0.0,1) 
Y=Y--YSP 

CALL  NUHBR(X,Y.YSCO, 0.0,1) 
Y=Y-YSP 

CALL  NUMBR(X.Y,ZSCQ, 0.0,1) 
X2=Xi+.25*XVMAX 
Yi=YVHAX 

CALL  BOXINtXi ,Y0 ,X2 ,Y1 , 0 .0,0.0) 
X0R=X2+XMARG 
YOR=YOR-fl .  CKYVMAX 
YSP=0.04*YVrtAX 
X  s  XQR 
Y  =  YOR 

CALL  LETTR<X,Y,4HY  =  ,0.0,4) 
Y=Y-YSP 

CALL  LETTRL  X,Y,4HX  =  ,0.0,4  ) 
Y=Y-YSP 

CALL  LETTR(X,Y,4HZ  =  ,0.0,4) 


Y=Y-1.S*YSP 

CALL  LETTR (X,Y,21H  PLT3P,  Ver.  3D:04:84  .0.0,21) 

X»X+DX3 

Y=YOR 

CALL  LETTR ( X , Y , HEDY ( 1 ) , 0 . 0 , NHEDY ) 

Y=Y-YSP 

CALL  LETTR (  X,Y,HEDX<1) ,0 . 0 ,NHEDX) 

YsY-YSP 

CALL  LETTR (  X,Y,HEDZ(i) ,0  0,NHEDZ) 

X0=X2 

X3=XVMAX 

Yi=YVMAX 

CALL  BOXIN(X2,Y0,X3,Yi,0 .0,0 .0) 

RETURN 

END 

SUBROUTINE  DATEP(X0,Y0) 

DIMENSION  IDATE(IS) 

CALL  FTIME(IDATE) 

CALL  LETTR <X0.Y0,IDATE<6), 0.0,20) 

RETURN 

END 

SUBROUTINE  REFAX ( XC , YC , ZC , XMAX , YMAX , ZMAX ) 

COMMON/DATA1 /XMARG ,XWMAX , XINF , YINF , XLlNF , YWMAX , XYMAR , ASPR 
COMMON/DAT A2/IPLX , IPLZ , THETA , PHI , DELTA , XSC , YSC ,ZSC , SIZE 
DATA  RD/S7. 295779513082/ 

DATA  HT0/.0i/,WTH0/0  001/ 

THETAR=THETA/RD 

PHIR=PHI/RD 

HT=HT0 

UHT=NHT0 

SINPH  =  SIN(  PHIR  ) 

COSPH  =  COS(  PHIR  ) 

SINTH  =  SIN (  THETAR  ) 

COSTH  =  COS(  THETAR  ) 

XLMX=XMAX 

YLMX=YMAX 

ZLMX=ZMAX 

SX=XLHX*XSC/ ABS ( XSC) 

SY=YLMX*YSC/ABS(YSC) 

SZ=ZLMXtZSC/ ABS ( ZSC ) 

XAXX  =  SX  *  COSPH 
YAXX  =  SX  4  SINPH  *  SINTH 
XAXY  =  0 . 

YAXY  *  SY  *  COSTH 

XAXZ  =  5Z  *  SINPH 

YAXZ  =  -SZ  t  COSPH  *  SINTH 

X  =  XC 

Y  =  YC 

CALL  ZMOVEI  X  ,  Y  ) 

X  *  X  ♦  XAXX 

Y  =  Y  +  YAXX 

CALL  ZDRAWI  X  ,  Y  ) 

X=X*HT 


CALL  LETTR(X,Y,iHX,0 .0  ,i) 

X  =  XC 

Y  =  YC 

CALL  ZMOYE<  X  ,  Y  ) 

X  =  X  +  XAXY 

Y  =  Y  +  YAXY 

CALL  ZDRAWt  X  ,  Y  ) 

X~X+HT 

Y=Y-HT 

CALL  L  ETTR ( X  >  Y , 1HY ,0.6,1) 

X  =  XC 

Y  =  YC 

CALL  ZMW£<  X  ,  Y  ) 

X  =  X  +  XAXZ 

Y  -  Y  +  YAXZ 

CALL  I DR AW (  X  ,  Y  ) 

X=X+HT 

Y=Y-2.*HT 

CALL  LETTR (X,Y,1HZ, 0.0,1) 

X  =  XC 

Y  =  YC 

CALL  ZMOVE<  X  ,  Y  ) 

RETURN 

END 

SUBROUTINE  BOXIN ( XC , Y 0 , XMAX , Y MAX , XYMAR . ASP R ) 

XrtAR=ASPR*XYMAR 
YMAR=XYMAR 
Xl=Xfl+XMAR 
Yi=Y0+YHAR 
X2=XMAX-X«AR 
Y2=YMAX-YHAR 
CALL  ZNOVEtXi ,Yi ) 

CALL  ZDRAW'Xi ,Y2) 

CALL  ZDRAW(X2,Y2) 

CALL  ZORAW ( X2 , Y i  ) 

CALL  ZDRAW(Xi,Yi  ) 

RETURN 

END 

SUBROUTINE  TITLE (XTMIN,XTMAX,YTHIN,YTMAX) 

COMMON/DATAi /XNARG , XWMAX , XINF , f INF , XlINF , YWMAX , XYMAR  t ASPR 
C0K«QN/DATA2/IPlX,IPLZ, THETA, PHI, DELTA , XSC , YSC , ZSC , SIZE 
CONMON/DATA4/ XMIN , XMAX , ZMIN . ZNAX , YMIN , YNAX 

COMNON/DATA9/NTITLi ,TITLEi(2) ,NTITL2,TITLE2(i5) ,NTITL3,TITLE3(iS) 
DATA  AN0/2HN0/AYE/2HYE/ 

YSP=0.5*YTMAX 

XG-XTHIN 

Y0=YTHIN 

XOR=XO+XMARG 

YOR=YTMAX-XMARG 

X=XOR 

Y=Y0R 

CALL  LETTR (X , Y ,TITLE2( i ) ,  C .  0 ,NTITL2) 

Y=Y-YSP 

CALL  LETTR <X,Y,TITLE3(1),0.O,NTITL3) 


Y=Y— YSP 

CAL L  BOXIN(XO , YO ,XTHAX ,YTHAX,0 . 0,0.0) 

RETURN 

END 

SUBROUTINE  DEVON 

C0MM0N/DATA6/INSTN , I FAST . IPLQT , IPMENU , IPEN<  2  > 

DIMENSION  SUBR ( 3 ) 

DATA  SUBR/4H  ,4H  ,4H  / 

LUT=i 

LUP=35 

90  WRITE ( 1 > 1 095 ) 

1095  FORNATrmmmtmmmmmtmmmtmmmtmmtm' 

C  ,/,'«*»  SPECIFY  THE  PEN  COLOR  FOR  THE  PLOTTER  **** 

C  ,/,"««  i.  BLACK,  2.  BLUE,  3.  RED,  4.  GREEN  *»**' 

C  ,/,"****  PEN  COLOR  =  ’) 

READ(1,»)  IPEN(2) 

IF  < IPEN < 2 )  LE  0)  IPEN(2)=1 
URTTtd.liOO) 

lioo  format; 

C  ./,“*«*  SPECIFY  THE  PEN  SPEED  FOR  THE  PLOTTER,  (1-35)  «#* 

C  ,/,'****  RECOMMENDED  VALUE  FOR  QUALITY  PLOTS  =  5  *«*' 

C  ,/.“*<**  PEN  SPEED  =  *) 

READd,*)  IPENd) 

IF  dPEN(l)  LE.  0)  IPENd >=5 
CALI  7PECN 

CALL  7AINT(LUT. IERROR) 

I?  i  IERROR. EQ.l)  GO  TO  100 
CAi.L  ERRORS(LUT . IERROR , 6HZAINT  ) 

GO  TO  999 

ICC  CALL  ZKINTdUT. IERROR) 

IF (IERROR  ECU)  GO  TO  200 
CALL  ERRORS(LUT, IERROR .6HZKINT  ) 

GO  TO  999 

200  CAu.  ZBINKLUT, IERROR) 

IF  (IERROR  EC*  3)  GO  TO  400 
CALL  ERRORS(LUT,IEAROR .6HZBINT  ) 

GO  TO  999 
400  ICQTP=C 

CALL  ZDINT(LUP,ICOTR, IERROR) 

IF ( IERROR. EQ.0)  GO  TO  600 
CALL  ERRORSCLUT, IERROR, 6HZDINT  ) 

GO  TO  999 
600  WRITE(1 ,101) 

101  FORMAT  ("  Type  <  /  )  to  continue  when  ready  =  "> 

READ  (l.*)  XN 

CALL  ZOESC(2050 ,2,0, IPEN , RDUM , IERROR ) 

IF  (1ERRCR  EQ.0)  GO  TO  999 

CALL  ERRORS ( LUT , IERRCR , 1 OHZOtSC  2050) 

999  RETURN 
END 

SUBROUTINE  ERRORS (LU , IERROR , SUBR) 

DIMENSION  SUBR (3) 

CALL  ZMCUR 

WRITE ( LU ,100)  IERROR, SUBR 


iflO  FORMATC<<<<<  Error  * ,15 , “  occured  in  subroutine,,)3A4,*>>>>>l,) 
RETURN 
END 

SUBROUTINE  VIEWP  (XUNAX YHMAX  .SIZE ,  I  ERROR ) 

DIMENSION  AR(2) ,VIEU(4) 

CALL  ZIUS(2S4, 9 ,2,IDUM, AR , IERROR ) 

XSIZE=i . 8 
YSIZE-AR<2) 

CALL  ZASPK<XSIZE,YSIZE) 

UXMIN*0 . 0 
UXMAX=XWMAX 
UYMINsO.OO 
UYNAX=YUNAX 

CALL  ZUIND <UXMIN , UXMAX ,WYMIN, WYMAX ) 

CALL  Z I WS (451 , 0 , 4, IDUM, VIEW, IERROR ) 

UNINXa . QStUlEU  <  2 ) /SIZE 
VNAXX=.954VIEW(2)/SIZ£ 

VMINY=  QStVIEM<4)/SIZ£ 

VMAXY  = . 9S4VIEU  <  4 ) /SI ZE 
CALL  ZVIEU  <  VMINX , VNAXX , VNINY , VMAXY ) 
XCSIZ=0.Qi5*(WXMAX-UXMIN)/SIZE 
YCSIZ=0 02S*<WYMAX-WYHIN>/SIZE 
CALL  ZCSIZ(XCSIZ,YCSI2) 

RETURN 

END 

SUBROUTINE  VPNAX(VIEW) 

DIMENSION  AR<2),VIEU<4> 

CALL  ZIWS <  451 ,0,4,1 DUN , VIEW , IERROR  > 

CALL  ZIU5<2S4,Q,2, IDUM , AR , IERROR ) 

XNEU=i.O 
YNEW=i . 0 

IF (AR(2)  LE. i .0)  YNEW=AR<2) 

IF(AR(2).GT.i.O)  XNE«=i./AR(2) 

CALL  ZVIEW(0.O,XNEW,O.O,YNEW) 

RETURN 

END 

SUBROUTINE  DEVQF 
CALL  ZCOLR(O) 

CALL  ZAEND 
CALL  ZBEND 
CALL  ZKEND 
CALL  ZDEND 
CALL  ZLEND 
CALL  ZEND 
RETURN 
END 

SUBROUTINE  LETTR  <X,Y,IAR,ANG,NC) 

DIMENSION  IA8(i>,VIEW<4> 

CALL  ZMOVE(X,Y) 

CALL  VPMAX(VIEU) 

CALL  ZTEXT(NC,IAR) 

CALL  ZVIEtKVIEUU) ,VIEW<2> ,VIEW<3) ,VIEW<4) ) 

RETURN 

END 


46 


SUBROUT INE  NUHBR ( X , Y , FPN , THETA , IFMT  > 
DIMENSION  JFMT(5) ,ISYM(3) 

GO  TO  (10,20,30)  IFMT 
10  ENCODE ( 1 5 , 15 , 1 S YM >  FPN 
GO  TO  100 

20  ENCODE! 15 , 25 , ISYH )  FPN 
CO  TO  100 

30  ENCODE (10, 35, ISYM)  FPN 
100  CALL  LETTR(X,Y, ISYM, THETA, 15) 

15  FORMAT(FS.l) 


25  FORMAT (1PE12. 4) 
35  FORMAT (F5.0) 
RETURN 
END 
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